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INFORMATIQUE 


INTRODUCTION 


Ce livre couvre l'étude des principales fonctions de l'électronique. Les 
programmes présentés ont tous été longuement testés et vérifiés. Cer- 
tains d’entre eux ont un aspect pédagogique ; d’autres permettent 
l'étude systématique de phénomènes fondamentaux : régimes transi- 
toires de circuits, réponses de filtre. Enfin, dernier aspect, le programme 
sur les boucles à verrouillage de phase constitue une introduction à la 
simulation numérique de systèmes non linéaires : étude passionnante, 
éliminant la complexité mathématique du problème. Nous pensons que 
le lecteur pourra aisément modifier tous les programmes proposés pour 
les adapter à ses besoins et à d’autres phénomènes. Les algorithmes 
employés, souvent très simples, sont toujours décrits avec précision. 
Ils ne doivent pas être considérés comme des ‘’modèles”’ parfaits, mais 
plus simplement comme des guides de travail. Dans un souci pédago- 
gique, ils sont toujours exploités de façon très claire. 


Le livre se divise en sept parties. La première partie aborde les filtres 
passifs ayant pour dénominateur commun le filtre LC ‘’bouchon'” ou 
résonnant (du 2° ordre) : circuits décalés, circuits couplés,.… circuits 
pour lesquels une étude mathématique directe est d’une complexité 
désarmante. 

La seconde partie s'intéresse au théorème de Fourier et à ses applica- 
tions. || est d’abord proposé une visualisation des formes d'onde de 
signaux électroniques classiques obtenues à l’aide d’un nombre de plus 
en plus grand d’harmoniques. Les applications du théorème sont ensuite 
exposées. Enfin, une méthode rapide de calcul des coefficients d’une 
série de Fourier est proposée sur un exemple. La troisième partie pro- 
pose l'étude des régimes transitoires des circuits les plus courants du 
premier et deuxième ordre. 

La quatrième partie aborde le filtrage ; filtrage du deuxième ordre, très 
classique, puis filtrage d'ordre plus élevé (jusqu'au 5° ordre). 

La cinquième partie est consacrée, dans un premier temps, à la recher- 
che d’un ‘’modèle’’ de la diode à jonction, tenant compte de certaines 
de ses imperfections. Quelques circuits à diode sont ensuite décrits 
(détection, redressement). Dans un second temps (détection 
synchrone), il est montré que trois dessins valent mieux qu'une longue 
formule. 


La sixième partie est en fait une longue introduction à l'étude des bou- 
cles à verrouillage de phase, composants essentiels dans les problè- 
mes de communication moderne, mais d’une complexité d'étude sans 
égal. Il est montré qu'une simulation numérique résoud bien des 
problèmes. 


Chapitre I 


RÉSONANCE 


CIRCUIT RÉSONNANT LC PARALLÈLE 


Ce circuit est extrêmement répandu en électronique. Il est à la base 
de la sélectivité des récepteurs radio. Son étude mathématique est sim- 
ple. Nous nous intéressons aux courbes donnant le module et la phase 
de sa transmittance. Les définitions à ce circuit sont les suivantes : 





On à : Z=—— avec 
1+jQ{(x — 1/x) x=F/Fo 


On a aussi : $= — Arc tg Q(x— 1/x). 


Nous proposons de tracer les deux courbes relatives à ces grandeurs 
en fonction de x. Nous normaliserons Z (en le divisant par R). 
L'échelle sur l’axe des fréquences est linéaire : nous la graduerons de 
x=0,5 à 1,5. Il est possible d'étudier plus finement ces fonctions en 
réduisant cet intervalle. 


LE PROGRAMME 

Lignes 60 à 120 : initialisation des variables. 

Lignes 130 à 180 : tracé des courbes sur l’écran. 
Sous-programme 210 à 290 : tracé des axes. 

Lignes 300 à 350 : choix de l’agrandissement des courbes. 


LES VARIABLES 


Q : coefficient de qualité du circuit étudié (ici, Q = 20). 

x _: Variable réduite de la fréquence ; x=f/fo. Dans cet exemple, x 
varie de 0,5 à 1,5. 

: module de la transmittance. 

: phase de la transmittance. 

: variable de la position horizontale. 

: facteur de l'échelle, valant successivement 1, 2, 4, 8, 


>CUN 


DÉROULEMENT DU PROGRAMME 


L'utilisateur doit d’abord donner la valeur du coefficient de qualité Q 
de son circuit. Les courbes donnant le moule (en haut de l’écran) et 
de l'argument (en bas) sont ensuite tracées. L'intervalle étudié va de 
0,5 fo à 1,5 fo. Deux petits traits, sur l’axe horizontal, rappellent les 
positions théoriques des extrémités de la bande passante ; ils sont donc 
situés à +fo/2Q de l'axe vertical. La valeur maximale de (Z) est 1 ; 
la phase varie de +7/2 à —x/2. 

Un INPUT permet de poursuivre le programme. Il est maintenant possi- 
ble d'agrandir les courbes par un effet de loupe. 

En répondant ‘‘O'’ à la question posée, un nouveau tracé de courbes 
va se produire, mais avec un agrandissement égal à un facteur 2. On 
obtiendra ainsi successivement 2, puis 4, puis 8 comme facteurs 
d'échelle. 


AGRANDISSEMENT : 


VOULEZ VOUS AGRANDIR LES COURBES? O 


AGRANDISSEMENT: 2e 


VOULEZ UOUS AGRANDIR LES COURBES? O 


10 REM HN 
20 REM +# 4 
JD REM *# CIRCUIT RESONANT # 
40 REM *# # 

50 REM HR HIERRRIERREREE 
70 CLS 

80 INPUT "Valeur de Q=";û 
90 A=! 

100 CLS 

110 U=iÙ 

120 GOSUE 210 

130 X1=1-0.5/A:X2=1+0.5/A 
140 FOR X=X1 TO X2 STEP O.C05/A 
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150 2=1/(1+(0æ(X-1/2))42)40.5 
160 P=ATNIQK(X-1/X)) 

170 PLOT 2#U,200+20047:PLÔT 244, 100+504P 

180 U=U+1:NEXT X 

200 GOTO 300 

210 REM AXES 

220 PLOT 20,200:DRAUR 420,0 

230 PLOT 20,100:DRAUR 420,0 

240 PLOT 220, 10:DRAUR 0,390 

250 PLOT 220+(A#200/0),195:DRAUR 0, 10 

260 PLOT 220-(4#200/0),195:DRAUR 0,10 

270 PLOT 220-(A#200/0),95:DRAUR 0, 10 

280 PLÔT 220+(A#200/0),95:DRAUR 0,10 

290 RETURN 

300 REM AGRANDISSEMENT 

310 LOCATE 20,1:PRINT "VOULEZ VOUS AGRANDIR®:LOCATE 20,2: INPUT LES CO 
URRES (0/N)°:4$ 

320 1F A$="0" THEN A=A#2 ELGE END 

330 GOTO 100 8 





CIRCUITS DÉCALÉS 


Le circuit résonnant LC a un accord pointu. Il est intéressant pour sélec- 
tionner une fréquence parmi d’autres. Pour certaines applications, de 
type passe bande, il a une courbe de réponse trop aiguë et les flancs 
pas assez raides. |l est alors posible d'utiliser deux circuits résonnants 
en cascade, à fréquences d'accord différentes, pour obtenir ce genre 
de gabarits. 

Dans ce cas, la transmittance globale au produit des transmittances 
correspondant à chacun des deux circuits (dont nous supposons les 
coefficients de qualité Q identiques). 





accord F1 7» accord F2 





En appelant X1=f/f1 et X2=f/f2 avec f1 et f2 pour fréquences de 
résonance. 

Nous poserons, pour simplifier : 

X1=X(1+N) et X2=X(1 —-N) ; X reste donc la variable réduite de la 
fréquence. Nous décalerons les deux circuits concernés de part et d’au- 
tre de fo, leurs fréquences de résonance respectives étant alors : 
F1=Fo/(1+N) # #Fo(1 —N) et F2=Fo/(1-N) # # Fo(1 +N) car N est 
petit. 

En faisant varier N, nous ‘‘’écartons’’ les deux courbes de façon symé- 
trique par rapport à Fo. 

Des programmes complémentaires montrent d’une part, comment se 
construit T à partir de T1 et de T2 et d'autre part, comment varie l’ar- 
gument de T (de +180° à — 180°) en fonction de la fréquence. On 
pourrait en déduire la variation du temps de retard de groupe du 
montage. 


LE PROGRAMME 


Lignes 70 à 110 :initialisation du problème. 
Lignes 120 à 240 : tracé des courbes correspondant à diverses 
valeurs de N ; O ; 0,005 ; 0,010 ; … ; 0,025. 
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Lignes 250 à 310 : tracé de la courbe correspondant à un circuit 

résonnant seul de mêmes caractéristiques : Fo, Q. 
Lignes 340 à 380 : possibilité d'agrandissement des courbes. 
Lignes 390 à 440 : dessin des axes. 


LES VARIABLES 


N : variable correspondant à l’écart entre les fréquences d'ac- 
cord des deux circuits. Cet écart vaut 2Nfo. 

Q : facteur de qualité des deux circuits. 

X, X1, X2 : variables réduites de la fréquence (X=f/fo ; X1=f/f1 ; 
X2=f/f2). 

Tiet T2 : modules des transmittances de chacun des deux circuits. 

T : module de la transmittance obtenue en mettant en cas- 
cade les deux circuits résonnants décalés. On a évidem- 
ment T=T1 T2. 

U : variable de position horizontale sur l'écran. 

P, P1, P2 : déphasages introduits par chacun des circuits concernés. 


A : facteur d’agrandissement des courbes. 


DÉROULEMENT DU PROGRAMME 


Il est identique à celui du programme : CIRCUIT RESONNANT LC 
PARALLELE. 

Les 3 variantes (programmes 1 à 3) permettent d'étudier les courbes 
de réponse (1), la phase (2), la construction des courbes de réponse 
(3). 





Intervalle étudié = 0,2 Fo autour de fo. 





Courbes obtenues avec le programme 1 : 


1 : circuit LC simple (fo) 

2 : 2 circuits LC identiques en cascade (fo) 

8 : 2 circuits LC décalés de 0,01 fo 

4, 5, 6 : 2 circuits LC décalés de 0,02 fo ; 0,03 fo ; 0,04 fo ; O,065 fo 














PROGRAMME 2 
N=0,03 





37 








#" 




















__ 1 É RS _ 


On visualise simultanément la courbe de réponse d’un filtre du deuxième 
ordre (courbe 2), celle de circuits décalés (figure 3) et la courbe de phase 
(courbe 1). 





PROGRAMME 3 
Construction de la courbe de réponse de circuits décalés : on visualise 


simultanément les courbes de réponse des deux circuits puis séparé- 
ment et la courbe résultant de leur mise en cascade. 


PROGRAMME 1 





ÎC RM RRéRRRFER#RERERERESE 
EG REM # x 
JC REM # CIRCUITS DECALES 4 
40 REN # # 

S0 SEM #RRHRERERESEHEREEREEE 
75 ÊLS 

BG INPUT “Valeur de G=":û 
90 &=! 

io6 CLS 

150 SOCUE JSC 

LéC FOR N=S Tü 9,025 STI QC 
136 U=û 

140 XA=1-0,1/A:X2= 1 +0. 1/4 


_ 
an 













4e 
2 
‘4e 
ae ‘dc 
2 
2 
< 
< 
# 
N 1e 
3 
£ = * 
Se: > du 
æ en eu La 
—. ll ci > # 
pra De ut Li 
SL <T < La 
«Z mnt Œ # 
Li OO  Æ 
in Œ Li Li 
Lu O LE 
Lo] Fi QU ci Œ 
ui Lx > LE Us LS 
C4 = F1 + = Fu 
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150 FOR X=XA TG X2 STEP G.00!/A 
LEO XI=K#11+N) 

165 Pi=-ATNIQEUXI-1/X1 5) 

170 T1=1/{{+{QætX1-1/411)A2)A0.5 
180 X2=X#(1-N) 

185 P2=-ATN(QK(XZ-1/#2)) 

196 T2=1/(1+(Qx(X2-1/Y2))/23 4.5 
ZCG T=TI4TZ 

Zi 
214 P=(PI+P2I4365 

215 PLOÔT PAU CP 

220 V=U+{ 

230 NEXT À 

250 U=G 

270 FOR X=XA T0 XE STEP O.061/A 

28 T=i/(i+(Qe(X-1/H1)425/005 

290 PLOT 20+2#Ù,10+380e7T 

305 U=u+i 

310 NEXT X 

340 REM AGRANDISSEMENT 

JC LOCATE ZO, {PRINT "UQULZZ VOUS AGRANDIR" :LOCATE 20,2:INPUT LES Cû 
URBES (G/N)° AS 

370 IF A$="0" THEN A=A+*2 ELSE END 

380 G0T0 190 

390 REM AXES 

409 PLOT 10,10:DRAWR 430,0 

410 PLOT Z20+ (4450 C/9},10:DRAWR 0,9 

420 PLOT Z20-(A#500/0),10:DR4UR 0,9 

430 PLOT 220, 10:DRAUR 0,395 

440 RETURN 





PROGRAMME 3 


LC REM #4 RO RER ES 
20 REM + + 
30 REM # + 6 RCUITS DECALES # 
40 REM * # 

SO REM RRRARRRREREEERREEE 
75 CLS 

80 INPUT "Valeur de ÿ="i 
90 Az! 

100 CLS 

{i0 GOSUE 590 


140 XA=1-0.1/4:XE=1+0.1/À 

150 FOR X=XA TO XB STEF C.GG1/A 

16G X1=XK CIN: 

170 Ti=1/(1+iqe(yi-i/X15HAZIAG,S 
180 X2=X#{1-N) 

190 T2=1/(1+(QKiX2-1/X2)342540,5 
200 T=TIKTZ 

210 PLOT 20+24U, 1C+280#7T 

211 PLOT 204240, 10+380#71 

242 FLOT 204240, 10+980#12 

220 U=U+] 

230 NEXT X 

340 REM AGRANDISGSEMEN 

360 LOCATE 20,1:PRINT “VOULEZ VOUS AGRANDIR" :LOCATE 20,2: INPUT ŸLES CO 
URBES (G/N)":A$ 

370 IF A$="0" THEN A=AK2 ELGE END 
380 G0TG 100 

390 REM AXES 

400 PLOT 10,10:DRAUR 430,0 

410 PLOT 220+(A4500/0),IG:DRAWR 0,7 
420 PLOT 220-(4#500/0), 10:DRAUR 0,9 
430 PLOT 220, 10:DRAUR 0,37€ 

440 RETURK 


CIRCUITS COUPLÉS 


La transmittance de circuits couplés est difficile à établir et délicate à 
interprêter. Il est cependant possible d'obtenir les courbes de réponses 
correspondantes par la méthode d'Evans. Cette méthode fait interve- 
nir les zéros-valeurs de jw qui annulent le numérateur de la transmit- 
tance (qui est une fonction complexe de jw) et les pôles valeurs de jw 
qui annulent le dénominateur. Ces valeurs particulières de jw étant trou- 
vées, tracer la courbe de réponse du système, revient à exploiter des 
considérations géométriques dans le plan complexe. Dans cette étude, 
nous ne nous intéressons pas à la phase de la transmittance. 

La méthode proposée déborde largement du cadre de l'étude des cir- 
cuits couplés. Elle s'applique à tout système dont on connaît les pôles 
et les zéros. Elle est très utilisée en servomécanisme. 

Notre étude portera sur le circuit suivant : 





Dans ce cas, on démontre que : 


D Me RO ST 
TT ac+Cn * (ppt) (p-p2 


Relation dans laquelle : 


CT 


k est le facteur de couplage =—— 
C+CT 


w0 est la fréquence de résonance des circuits R, L, C pris séparément 
et supposés ici identiques. Q est leur facteur de qualité commun. 


p1 et p2 sont les pôles de la fonction : on a : 


wO kKwO . kwO 
4 he PRO URT Se EN nice 
pi 2Q ! J 2 ie 20 ju0— ] 2 
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p = variable complexe = jw. 


Dans le plan complexe, on a : 


— wo / 2Q 


PARTIE 
IMAGINAIRE 


re) Iwo + | kwo 
2 


PARTIE 
RÉELLE 





On constate que le module de Zt peut s'écrire : 


kKwO 1 


Fr. 4(C + Ct) ' r1 r2 

en appelant r1 et r2 les distances séparant les pôles du point jw. r1 et 
r2 sont les hypoténuses de triangles rectangles faciles à exploiter. jw 
correspond à la pulsation pour laquelle on calcule Zt. 

Pour avoir l'allure des courbes en fonction de k, on peut se contenter 
d'évaluer k/r1 r2. 

Nous étudierons un circuit pour lequel Q = 10, fo=1 MHz. k prendra 
des valeurs comprises entre 0,01 et 0,8. La fréquence f variera de 
0,5 MHz à 1,5 MHz. 


LES PROGRAMMES 


Lignes 70 à 110  : entrée des données. 

Lignes 150 à 240 : tracé de la courbe. 

Lignes 270 à 350 : sous-programme de tracé des axes et graduations. 
Lignes 360 à 390 : possibilité de ‘‘loupe’’ sur les courbes. 


Le programme n° 2 diffère du programme n° 1 par le fait que les cour- 
bes correspondant à plusieurs valeurs de k prédéfinies (ligne de DATA) 
sont affichées simultanément. 
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LES VARIABLES 


fet Q : tels que définis plus haut. 

f : fréquence considérée. 

k : facteur de couplage (son nom est G en variable de 
boucle). 

A et YO : facteurs d'échelle. 

R1 et R2 : distances correspondant aux deux pôles. 

B, C, D, E, U, N : noms donnés à des constantes. 

F1, F2 : limites de la bande de fréquence étudiée. 

S : pas d’incrémentation sur l’axe des fréquences (il y 


a 200 points explorés). 


DÉROULEMENT DU PROGRAMME 


Celui-ci affiche les courbes de réponse de circuits couplés, soit pour 
une valeur donnée de k (programme 1), soit pour une suite préétablie 
de valeurs de ce paramètre (programme 2). 

L'axe horizontal est gradué. Lors du premier passage, le facteur d’agran- 
dissement A vaut 1. La zone de fréquences explorée a alors une lar- 
geur de 10 FO/Q, centrée sur FO. Le pas d'incrémentation S vaut alors 
1/200 de cette valeur. Il y a 10 graduations sur l'axe horizontal, espa- 
cées donc de FO/Q. 





Légende de la photo (programme 2) ; Q=10. 














PROGRAMME 2 
(k varie de 0,01 à O,8). 


Chaque passage par une réponse positive à la question ‘voulez-vous 
agrandir les courbes ?'’ se traduit par une dilatation d’un facteur 2 de 
l'échelle horizontale par rapport à celle qui lui est immédiatement pré- 
cédente. 





PROGRAMME 1:0= 10 k=0,2 


PROGRAMME 1 

10 RENM ee RE 
20 REM # 5 
JC REM # CIACUITS COUPLES +# 
40 REF * F 
SO REM RARE ERREUR 
70 €LS 

86 Sc 





150 GÜCUE 270 
1éê R- Ne ones 
170 Fi=FC4{i-15/(GxA))) 
180 F2=F0#(1+(5/{5#4))) 
190 G=16+F0/{2C0#4#0) 
200 FOR F=Fi TO F2 STEP S 
210 RI =ÜCARS(F-FCHFOR TC: 

220 RZ ={(ABS(F-FO-FÜ4K/D))AD+(E 
230 N=N+{:PLOT CN, LOHZ#(C,G/CRIXE 
2460 NEXT F 

250 GÜTO 360 
270 REM AXES 
272 LOCATE JO, {:PRINT 'Q=*:0 
273 LOCATE eo "Kat :K 
274 LOCATE 10,5:PRINT ‘Foz": 
280 PLOT Z0,20:DRAUR 400,0 
290 FOR X=220 TO 420 STEE 4G+A 
JO0 FLOT X,20:DRA#R 0,12 
310 NEXT 
320 FOR X=220 TO 10 STEP -40#A 










k. 


Le 


cm 


JC PLOT X.2C:DRAUR 0,10 


34G NEXT 

350 RETURK 

360 REM AGRANDI: 
370 LOCATE 1,2841N É 
JeC IF A$="0" OR A$=* 
90 GÜTO 130 








OURBES (ON: AS 





Ce) 
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PROGRAMME 2 


IC REM RERANÉREEERRRÉEERRRE 
Z0 REM +# * 





Ce 
Det] 


Ces 
“0 
on 
Lg 
+ 
# 


REM HEERIEREREFERRRREREESE 


a 


LLS 


CO —J IN + 
Let] 


C2 


-Q 


1 0 CD -41 Ce Un Es tn ve CD C2 











ge are 
Ces] 


C2 CI 


em | 


Où 1 


Nue ne pe et pue en 


290 N=N+ 
240 NEXT 


D4N,2042# (0 


280 PLOT 20, 20:DRAUR 400,0 
290 FOR X= 


L 


TO 420 STEP 4044 








J00 FLÔT X,20:DRâNR Q, it 

315 NEXT 

330 EF -à0#é 
330 PLOT K,ZC:DRAMR 9.10 

340 NEXT 

355 RETURA 

d60 REM asso [T 





; AT a+ 
LÜCAT E trs 


HEGEMENT LES COUREES (GNT 8E 


Ci 
a 






N Azhx: ELSE EXD 





Æ Li Ci 
C3 “0 CCI +4 
EN C9 C1 

re 

— 


Chapitre II 


FILTRAGE 


FILTRES PASSIFS 


Nous nous limiterons au cas des filtres passifs symétriques non dissi- 
patifs, fermés sur leur impédance caractéristique (Zc). Leur schéma 
général est le suivant : 


FILTRE ENT FILTRE EN TT 





La transmittance d’un filtre passif est étudiée en deux étapes. D'abord, 
calcul d'une courbe universelle, valable pour tous ces filtres, ensuite 
l'application au cas particulier : ici, un filtre passe bas, dans lequel Za 
est une inductance et Zb un condensateur : 


On démontre que pour ces deux circuits : 


T = RE k+ T2 VR (K+2) ; —o<k< +00 


alors : 

si—-2<k<0:T=1 et Arg (T)=Arc Cos (1+k) 
sik<—2:T=(1+k)+Vk (k+2) ; Arg (T)=7 

sik>O:T=(1+k)-Vk (k+2) ; Arg (T)=0 

Nous prendrons comme exemple le filtre en T, dans lequel Za=jLw et 


Zb=1/jCw. 
Alors k= —LCuw? varie entre — © et O (on peut écrire : k= — (f/fo)?). 


LES PROGRAMMES 


Le programme n° 1 permet la construction de courbes universelles en 
fonction de k. Le programme n° 2 fait correspondre, pour un circuit 
donné, module et phase de la transmittance aux valeurs de f envisa- 
gées. 
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Programme n° 1 


Les valeurs du module T et de la phase P de la transmittance sont cal- 
culées par zones de valeurs de k : 


Lignes 90 à 140 


Lignes 140 à 210 : 


Lignes 220 à 260 
Lignes 280 à 430 
Lignes 440 à 450 


Programme n° 2 


: —10<k<-2 
—1,9<k<-0,1 
: O<k<10 


: tracé des axes et des graduations. 
: sous-programme d'affichage des points calculés. 


Il permet de tracer directement module et argument de la transmittance 
dans un cas particulier. Ici, nous avons : K = —(F/F0)? (voir le texte) ; 
la fonction est étudiée de F=Fo à F=2F0. 


LES VARIABLES 


Programme n° 1 


| ou X T7 


, NL 
Programme n° 2 


T (200) et P (200): 
F 
X 
N 


: paramètre utilisé dans les fonctions T et P. 

: module de la transmittance. 

: phase de la transmittance. 

: variable d’incrémentation, assurant la position hori- 


zontale des points sur l'écran. 


: Variables de boucles locales (tracé d’'axes.…. }). 


mémorisation des différentes valeurs de T et P. 


: fréquence considérée. 
: variable d'incrémentation. 
: Valeur de l'indice des matrices T(200) et P(200). 

















PROGRAMME 1 


10 SEM HER ÉSÉRÉRERTE 


20 REM * * 
30 REM # FILTRES PAGIFE # 
40 REM + L 


50 REM #KHRHERRRÉNÉRREREE 

60 INK CO: INK 1,Z6: PAPER CiPEN 
70 CLS:X=0 
BG GGSUR ZBO:REM AXES 

G FGR K= -10 TG -3 STEP 0,1 

100 T=ARS(K+i+(HAZ4C4K) AC, Si:B=j85 
110 GOSUE 440 


30 NEXT K 

140 FOR K=-1,7 To -9,1 STEP 0.1 

150 iF K=-1 THEN P=SQ:GûT0 190 

160 7=1:P=CI80/PI#ATNE (RATE) /(K#1 7); 
170 IF K£-1.05 THEN P=i89-P 

190 ÏF Ki-1,05 THEN P=ARS (PF: 

190 GOSUE 440 

200 X=X+i 

210 NEXT K 

220 FOR K=0 T0 10 STEP G.1 


230 T=itK-({K42+2#K) 40,5) :P=5 


240 GOSUE 440 
250 X=X4# 
280 NEXT K 
270 END 
280 SEM AXCS 
290 PLÔT 20, 10:DRAUR 400,0 
300 PLOT 20,20C:DRAUR 400,0 
310 FOR U=D TO 95 STEP 9 
320 FOR N=0 TO 200 STEP 10 
330 PLOT 20+24N, 10+2#U:DRAUR C, 10 
340 NEXT N:NEXT U 
350 PLÔT 220, 10:DRAUR 0, 180 
360 PLOT 220,200:DRAUR 0, 180 
370 FOR L=10 T0 :80 STEP 9 
380 PLOT 220,L:DRAMR 10,0 
390 NEXT L 
400 FOR L=200 TO 370 STEP 29 
410 PLÔT 220,L:DRAUR 10,0 
420 NEXT L 
422 PLOT 220,191: DRAUR -5,-S:PLOT 220, 191 :DRAUR 
425 PLOT 220,381:DRAMR -5,-S:PLOT 220,381 :DRAUR 
427 LOCATE 15,1:PRINT "PHASE" 
428 LOCATE 15,14:PRINT “1 MODULE" 
29 LOCATE 16,5:PRINT "180" :LOCATE 1.1 
420 RETURN 
440 PLOT 2042#X, 10418047 :PLOT 20+24X, 200+24P/3 
450 RETURN 


PROGRAMME 2 


10 REM HR AE 


20 REM * x 
30 REM + FILTRES PASSIFS # 
40 REM x * 


50 REM HENRI 

60 INK C,O:INK 1,26:PAPER O:PEN 1 
70 CLS:X=Û 

100 DIM T(200):DIM P(200} 

120 FOR K= -10 70 -2 STEP G.i 

190 T=ABS(K+1+(KA2+2#K)4G,5) :P=i80 
140 TIX)=T:P(X)=P 

150 X=X+i 

155 GOSUR 600 


3, 
5 


5 
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160 NEXT K 
170 FOR K=-1.9 TO -0.1 STEP 0.1 
180 IF K=-1 THEN P=90:G0T0 220 
190 T=1:P=(10C/PIIRATNC (KA24Z#K)/(K41)) 
200 IF K<-1.05 THEN P=180-P 

210 IF K3-1.05 THEN P=ARS IP: 

220 T(X)=T:P(X)=P 

230 X=X+1 

235 GOSUR 600 

240 NEXT K 

250 FOR K=0 TO 10 STEP 0.1 

260 T=14K-C(KA2424K) 40.5) :P=0 

270 TUX=T:P(X)=P 

280 X=X+1 

285 GOSUR 400 

290 NEXT K 

295 PLOT 20,10:DRAWR 400,0 

300 CLS 

310 REM AXES 

320 PLOT 20,10:DRAUR 400,0 

330 PLOT 20,200:DRA4R 400,0 

340 FOR U=0 TO 95 STEP 95 

350 FOR N=0 TO 200 STEP !C 

340 PLOT 20424, 10+2#U:DRAUR 0, :0 
370 NEXT N:NEXT U 

380 PLOT 220,10:DRAMR 0, 180 

390 PLOT 220,200:DRAUR 0, 180 

400 FOR L=10 TO 180 STEP 18 

410 PLOT 220,L:DRAUR 10,0 

420 NEXT L 





423 PLOT 220,191:DRAUR -5,-3:PLÔT 229,191 


425 PLOT 220,381:DRAUR -5,-S:PLÔT 220,381 
427 LOCATE 15,1:PRINT "PHAGE" 

428 LOCATE 13,16:PRINT *! MODULE" 
429 LOCATE 16,5:PRINT "IC" :LOCATE 1,1 
430 FOR L=200 To 370 STEP 20 

440 PLOT 220,L:DRAUR 10,0 

450 NEXT L 

510 REM EXEMPLES 

520 FO=1 

530 FOR F=F0 TO 24F0 STEP FO/200 

540 K=-(F/FO)42:N=104K+100 

550 T=(TCINTEN) )#TCINTUN#E)))#0.5 
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:DRAGR 5,5 
:DRAUR 5,-5 





360 P=(PCINTIN))+PBCINTIN+I 5) ix0,S 
570 PLÔT 20+F 4200, 10419047: PLGT GOHP450O, 200+287/3 
S8û NEXT F 
370 LOCATE 1,i 
60G LÜCATE 1,5: 
&iG LOCATE 1,15: 
&àC RETURN 


END 
PRINT “PATIENCE... J'en suis 
PRINT *Et je vais jusqu'a 2GC ‘l'* 








pui L. 


DUT AS ee ni 2 








PROGRAMME 2 


Etude du filtre par lequel k= —(F/F0}? ; F varie de Fo à 2Fo. 
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CONSTRUCTION D'UN FILTRE 
DU DEUXIEME ORDRE 


Il existe différentes possibilités de construire un filtre passe bas du 
deuxième ordre. On sait l'importance de tels filtres ; il est en effet pos- 
sible de réaliser un filtre du nième ordre par mise en cascade de filtre 
du second ordre (et éventuellement du premier ordre). Nous avons choisi 
de réaliser un filtre actif contenant une source commandée de tension. 
Il s'agit là d'un amplificateur de gain positif que l’on sait facilement réa- 
liser à l’aide d’un amplificateur opérationnel. 

Le circuit proposé est le suivant : 





La source commandée (K) peut être réalisée de la façon suivante : 





v' RER 


G : amplificateur opérationnel. On a alors Eva = Kz=1 
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R’ 


Pour le filtre proposé, on démontre que : 


V2 K 
Vi R1R2C1C2(ju)2 + [R2C2 +R1C2+(1—K)R1C1Jjw+ 1 


Un tel système est défini par : 

— son gain Go dans la bande passante (en dB) (= 20 log K) ; 
— son coefficient de surtension S (en dB) ; 

— sa bande passante F3 à —3 dB. 


Toutes ces quantités dépendent directement des valeurs de R1, R2, C1, 
C2. L'allure des courbes obtenues est la suivante : 


Gain en dB 


| S (en dB) 


Go 
G=sid—=— 
Lu OdB log (fréquence) 





Le programme propose donc de réaliser un tel filtre. 

Il'est d'abord demandé à l'utilisateur les valeurs souhaitées de Go, S, 
F3. À partir de ces données, le programme calcule des valeurs possi- 
bles pour C1 et C2, en s'étant imposé R1 =R2 = 10 KA. Il y a en effet, 
au départ, 4 inconnues (R1 R2 C1 C2) pour 3 données. On peut donc 
s'imposer R1 et R2. Il est ainsi proposé les valeurs de tous les éléments 
du circuit (K est imposé par Go). Si ce valeurs ne conviennent pas, l'uti- 
lisateur peut alors proposer ses propres valeurs de R1 et R2. Le pro- 
gramme en déduit de nouveau C1 et C2. 

Quand le choix des éléments est fait, le programme passe en haute 
résolution. Il dessine d'abord le gabarit correspondant au circuit désiré 
par l'utilisateur (Go, S, f3) et dessine la courbe théorique correspon- 
dante. 

Puis, il demande les valeurs effectivement employées pour les 
composants, sachant que ceux-ci sont pris ‘‘dans le tiroir’’, c'est-à-dire 
correspondant à des valeurs normalisées. Il dessine alors la courbe cor- 
respondante, ce qui permet à l'utilisateur de juger si le filtre obtenu lui 
convient, dans la mesure où seuls les points des deux courbes qui ne 
coïncident pas restent affichés. 
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Les axes représentés sur l’écran sont gradués en décades pour celui 
des fréquences, par bonds de 20 dB pour celui du module de la trans- 
mittance. Ainsi, l'axe OX part de Fo/100 et va jusqu'à 100 Fo (2 déca- 
des de part et d'autre de Fo) ; OY est gradué de +20 dB à — 80 dB. 
Les relations mathématiques utilisées sont les suivantes : partant de 
S, on en déduit l'amplitude de la résonance. 


S x 2,3 


TR= 
exp { 30 


) 


et de là, le coefficient d'amortissement du circuit : 


A) 1-V1-1/TR2 


2 


puis, 


= 2— V = 
U= NV TPE METRE qui vaut F3/F0 


d'où la valeur de la fréquence propre du circuit Fo =F3/U. 


| / Ym2=2{(1=Ki 
Enfin, posant de nouveau U = _ = D 


, 
et V=C1C2= ————— 
R1R2w0? 


on obtient C2=VU?V et C1= , Sachant que : 


(2V = +(1-KRV a 


On en déduit C1 et C2, connaissant K, R1, R2. 

Un sous-programme (630) recalcule les coefficients A et B de la fonc- 
tion de transfert à partir de R1, R2, C1, C2. Cette démarche est néces- 
saire puisque l'utilisateur peut faire intervenir ses propres valeurs de 
composant. 


m = 


Le 
2 


LE PROGRAMME 


Lignes 100 à 150 : entrée des données numériques. 

Lignes 160 à 250 : calcul des caractéristiques du filtre : m et FO. 
Lignes 280 à 320 : calcul de C1 et C2. 

Lignes 330 à 400 : affichage des valeurs conseillées pour les compo- 
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sants. 
Lignes 410 à 460 : modification possible des valeurs proposées par 
R1 et R2. 
Lignes 470 à 510 : tracés du gabarit, des axes, des graduations. 
Lignes 520 à 580 : tracé de la fonction (une décade occupant 50 
points, l’incrément est de 50010 = 1010,02 =D). 
Lignes 630 à 690 : introduction des valeurs numériques des éléments 
effectivement possédés par l'utilisateur, et dessin 
du tracé correspondant. 
730 : calcul des coefficients A et B. 
830 : sous-programme de tracé des axes et des 
graduations. 


Lignes 700 
Lignes 740 


D Q- 


LES VARIABLES 


Elles ont toutes été explicitées dans les paragraphes précédents. 


F3 : fréquence de coupure à —3 dB souhaitée. 
S : surtension souhaitée en dB. 
GO : gain souhaité en dB dans la bande passante. 
TRet2Z : correspondent à S et Go sur une échelle linéaire. 
V,U : variable temporaire de calcul. 
M : coefficient d'amortissement du filtre. 
Wo et Fo : pulsation et fréquence propre du filtre. 
R1, R2, C1, C2, K : valeurs des éléments du filtre. 
D : incrément pour les fréquences (= 5010). 
W: : valeur du module de la transmittance à la pulsa- 
tion W. 
A et B « coefficients de la fonction de transfert. 
V2 _ 1 


(— 


Vi 1 + Bjw + A(jw)? 





CRLCIIL D’'IUN FILTRE 
PASSE-BAS 
DU SECOND CRDRE 
Frequence de coupure x —-3daB :27 2E4 


Surtension souhaitee en dB<29>S>9):7 


" 
G « 


Gain dans la bande passante 
{ DS aB sonsei:le > : 


Introduction des données. 
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” 3 à 
et due dl A'amortissement du civrceu 
et sa frequence naturelle est 
G.13E+60SK7. 
Nous vous proposons les valeurs suiva 
S pour les composants: 
R1i= 1.60E+:04 OHMS 
R£= 1.G060E+04 OHMS 
C1= S.SBE+O8 nF 
C2= 2.44E-S2 KE 
Ii est ohiigatoire üe prenûäre K= 
Ces valeurs vous conviennent-ellses7 © 


Calcul des paramètres du filtre. 





ns LES ELEMENTS CONSEILLES? N 
La courbe de réponse théorique du filtre est tracée. 


Dans ce cas,quels sont les valeurs 


elements que vous possedez? 
R1i=7?7 1E4< 
R2=7 1E4 

1=2 6.8E-9 
C2=?2 9B.22E-9 
puis, suivant les composants dont on dispose. 
38 


des 





AUVEZ-UQOUS LES ELEMENTS CONSEILLES? © 


la courbe réellement obtenue est tracée. 


10 REM HER HOOREERENONER 
20 REM + 4 
3G REM # CALCUL D'UN FILTRE * 
40 REM # PASSE--HAS x 
50 REM * DU SECOND ORDRE  +# 


bD REM *# x 
70 REM HER OEM RES 
80 CLS 


120 LOCATE 10,2:PRINT "CALCUL D'UN FILTRE" :LOCATE 14,3: PRINT" PASSE-FAS 
iLOCATE 11,4:FRINT "DU SECOND ORDRE* 


130 LOCATE 1,6:INPUT ‘Frequence de coupure a -3dE 1373 

140 LOCATE 1,6:INPÜT “Surtension souhaitee en dB(9:0;0}:":6 

150 LOCATE 1,10:INPUT “Gain dans la bande passante { 0 dB 
conseille } :":G0 


160 TR=EXP(S#2, 3/20) 

170 U=(1-(1-(1/(TR42)))A0.5)/2 

180 M=UAD.5 

190 H=4xi-2 

200 U=(-H+((HA2)+4)40,5)40,5/(240.5) 
210 WO=24PI+FJ/UÙ 

220 FO=KkO/(2*81) 

230 K=60 

240 12=EXP(G0+2, 3/20) 

250 CLS 

260 LOCATE 1,É:PRINT * Le facteur d'amortissement du circuit est:”: 
USING °##, H#° M 


265 PRINT 

270 PRINT ‘et sa frequence natureile est E'USING FH HRAAAAR 
WO/(2KPI)::PRINT "Hz." 

280 R1=1C000:R2=16000 

290 U=(2+M+((4xMxM)+9#(17-15)40.5)/4 

300 V=1/(RI#R2# (4042) 

310 C2=(V#(U42))A5,5 

320 Ci=V/C2 

340 LOCATE 3,8:PRINT ‘Nous vous proposons les valeurs suivantes goùr | 
es composants!" 

330 LOCATE 10,12:PRINT "RI=' USING HKHHAAAAS RTS: PRINT * OHME® 

JG LÔCATE 10,14:PRINT "R2=":UGINS “#H,HAAAA SRG PRINT S CHMG* 

370 LOCATE 10,16:PRINT *CI="S;UGING MHHQHRAAAAR SCIE TACGS PRINT n° 
J80 LOCATE 10,{8:PRINT "CZ=" USING HA, HHMANAR CE IEACUS SPRINT # nF* 
390 PRINT:PRINT " Il est obligatoire de grendre K=":;USING *#H#.#":7 
400 PRINT:INPUT "Ces valeurs vous conviennent-elles":A$ 

420 IF A$="O* OR A$="o" THEN GGTO 470 

430 CLS:PRINT' Quelles valeurs proposez-vous?* 

440 PRINT: INPUT'RI ="3R1 

450 PRINT: INPUT'RZ =":R2 

460 GOTO 290 

470 CLS 

480 GOSUR 710 

490 WINDOUW #1,2,30,1,4-((5+K)/9):PABER #1,3 

495 FOR 1=2 TO JO:FOR J=i TG 5-(SrK)/16: PRINT #1," ":NEXT JiNEXT I 
500 AJ=(INT(95+50#L0G(F3/F0)/2,3))/8:B3=(48-K;/B:D3=(1504K) /8 

502 WINDOW #2,2,A9,B3,D3:PAPER #2,3 

505 FOR 1=2 TO AG:FOR J=B3 TO DA:PRINT #2, *:NEXT JiNEXT I 

510 GOSUR 740 

520 D=1040.02 

530 W=H0/(100*D) 

540 FGR X=X=0 TO 199 

550 Af=(1-AxW#4)A2:B1= (BW) AZ:W1=2/((A1+B1) A0, 
560 PLOT 2#X,322 + ((40#L0G(W1)1/2,3) 

570 W=WxD 

580 NEXT X 

630 LOCATE 1,24:INPUT "AVEZ-VOUS LES ELEMENTS CONSEILLES" ;A$ 

640 IF A$="C" OR A$="0* THEN END 

650 CLS:PRINT "Dans ce cas,quels sont les valeurs des elements que vo 
us possedez?* 

660 PRINT: INPUT'RI=":R1 

670 PRINT: INPUT*R2=" ;R2 

680 PRINT: INPUT'CI=":Ci 

« 


in 
— 


690 PRINT: INPUT" CZ= "3 CE 

695 GüTà 470 

706 REM CALCUL DES COEFFICIENTS 
710 A=RIXREHCI#CE 

720 B=iRIKCZ+RZ#CEHII-ThRI EC 
730 RETURN 

740 PLOT {,422:DRAKR 486,0 

750 LOT 1,20: DRAUR 0,385 

760 FOR N=i Tü 3 

710 PLOT 1CO*H,32G:DRANR 0,6 
790 NEÏT N 

800 FÔR N=û Tü 5 

BIC PLOT 1,405-40#N:DRA4R 8,0 
826 NEÂT N 

8230 RETURN 


a 


FILTRES ACTIFS DU 
DEUXIEME AU CINQUIEME ORDRE 


Un filtre est souvent caractérisé par la vitesse à laquelle chute le gain 
au-delà de la bande passante. Un filtre du 2° ordre coupe à raison de 
40 dB/décade ; cela veut dire que son coefficient de transmission est 
divisée par 100 chaque fois que la fréquence est multipliée par 10. Cela 
peut être insuffisant dans certaines applications ; un filtre d'ordre n 
coupe à 40 n dB/décade. De tels filtres sont en général réalisés par mise 
en cascade de filtres d'ordre 2 et d'ordre 1. 

Parmi l'infinité de réalisations possibles, il est classique de privilégier 
certains critères de choix ; l'existence de filtres normalisés s'ensuit. 
Nous avons choisi de présenter quatre de ces types de filtres (voir le 
programme). Le programme propose de dessiner les courbes de réponse 
de tels filtres et de donner l'aspect de la réponse en signaux carrés (voir 
les programmes sur FOURIER). 


LE PROGRAMME 


Lignes 80 à 220 : présentation des différents types de filtres. 
Lignes 220 à 300 : choix du type de filtre (Z) et de l’ordre (Y). 
Lignes 310 à 370 : tracé des axes sur l'écran. 

Ligne 380 : aiguillage vers le type de filtre. 


type Fe Dent nv H chebyscheffl  Bessel 
Z= Ru Ra 4 


D “410. 1 | 1350 | | 1820 | 


Lignes 2270 à 2320 : définition des fonctions. 

Lignes 2330 à 2430 : inscriptions sur l'écran. 

Lignes 2460 à 2530 : calcul des indices des matrices. A et B corres- 
pondant aux différents termes de la série de 
Fourier. 

Lignes 2560 à 2610 : inscriptions sur l'écran. 

Lignes 2650 à 2700 : calcul de la tension de sortie du filtre. 

Lignes 2710 à 2720 : calcul du temps de réponse. 

Lignes 2740 à 2760 : inscriptions sur l'écran. 
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ALGORITHMES MIS EN JEU 


TRANSMITTANCES DU PREMIER ORDRE 
_ 1 1 


Tai T}# 
1 +jw/w0 VT+ (w/w0)? 


Arg (T )=—Arc tg © 
wO 


TRANSMITTANCES DU DEUXIEME ORDRE 
1 





TT —  _— 
1+2jm = +(j =} 
wO wO 
(TE ; 
1 uw? 124 (_2M &0 12 
wO? &wO 
2Mw/w0 


A T )=-Arc tg——— 
d'ALE, FE DE ere 


(varie de O à —7) 


Cette dernière relation pose plusieurs problèmes lorsqu'elle est employée 
sur un micro-ordinateur. En effet, ceux-ci ne possèdent en général que 
cette seule fonction trigonométrique qui a l'inconvénient de n'être défi- 
nie qu'entre — x/2 et + x/2, et possède une disontinuité pour 7/2. 


Nous ajoutons ++/2 à Arg (T ) en écrivant : 


T'- 2jMmw/w0 


1+ 2jm = +} 
wO wO 


1 





pepe) 
2M  wO ra) 


1 & wO 
On a alors Arg (T')= — Arc tg. Ce { ee 2) 


mais —7n/2<Arg (T')<+7#+/2, ce qui ne pose plus de problèmes. 
On en déduit : Arg (T)= Arg (T')- x/2. 

Pour chaque cas étudié, la fréquence varie de F0/100 à 100 Fo. 
Partant de F0/100, chaque pas est obtenu en multipliant la fréquence 
par 50V10, soit 1/50° de décade. Au bout de 50 pas, on est ainsi à 
Fo/10,.. 

Pour chaque fréquence, le module de la transmittance est rangé dans 
A(200), la phase dans B(200). 

Un signal rectangulaire a pour développement en série de Fourier : 


n=o 1... ; ; 
y{t)= £, Sin 27 nFt avec n impair. 
n = 


Nous nous imposons comme fréquence de base Fo/4 (cette valeur sem- 
ble raisonnable). 

Pour l'harmonique d'ordre n, l'indice correspondant des tableaux A(200) 
et B(200) est donné par : C{n)= 50 log (25 n). 





[En effet, l'indice du tableau A (i) correspond à la fréquence 
f=(f0/100) —(10)0.02xi- 
: fo fo fo 
sif=n —,0on an— = x (10)0.02i 
4 4 100 (19) 


D'où 0,02i=log{(25 n) et donc i= 50 log (25 n) 
Le temps du retard 7 est donné par : 


Tr=d@/dw. 


La phase © est rangée dans B(200). Entre deux valeurs successives 
de F, on a donc : 


= = Bo) AF étant la variation de fréquence. 








27%AF 
On a donc LE x Log 10 x (10)9.02n x 0,02 avec dn=1. 
dn 100 

D'ais= = Bin+1)-Bi{n) 

27 x 0 x 2,3 x (10)0.02 x 0,02. 

100 

LES VARIABLES 
Z : type du filtre (1<Z<4) 
YetR : ordre du filtre (2<R<5) 


«4 


X : pulsation réduite. 


D : incrément = (10)10,02. 

A(200) : matrice contenant les 200 valeurs du module de la 
transmittance. 

B(200 : matrice contenant les 200 valeurs de l'argument de 
la transmittance. 

ENP(X) : fonction du 1°" ordre (module). 

FNC(X) : fonction du 1e" ordre (argument). 

ENSIX) : fonction du 2e ordre (module) 

FNQIX) * j 

NES : fonction du 2e ordre (argument). 

A, B,C, D,E,F : coefficients des fonctions. 

C(30) : matrice contenant les indices des matrices A et B 


correspondant aux différentes harmoniques de la 
série de Fourier. 


K : numéro de l'harmonique considérée. 
FNP(X) 1 

FNC(X) 1+jCX 

FNS(X) 1 

FNA(X) 1+)jBX + A(jX)? 

FNQ{(X) 1 

FNBIX) 1+jFX +E(jX)2 


Nous ajoutons + 7/2 Arg (T) en écrivant : 
REPONSE A UN ECHELON 


Dans ce cas, l'échelon est assimilé à son développement en série de 
Fourier limité à ses 15 premiers termes. Ceci est d’ailleurs justifié par 
le fait que les circuits étudiés sont du type passe-bas, et donc que l'in- 
fluence des harmoniques d'ordre élevé est très faible. La période de 
base du signal rectangulaire est 4/f0 ; le quinzième terme de la série 
Fourier a donc pour fréquence 31 fo/4 ; il est donc atténué au moins 
dans un rapport 100 (pour un circuit du 2e ordre) par rapport aux ter- 
mes situés dans la bande passante. 

Chaque terme de la série voit donc son module affecté de la valeur de. 
la transmittance du filtre étudié à la fréquence correspondante, et sa 
phase tournée de l'argument, d’où l'intérêt de calculer celui-ci pour cha- 
que fréquence lors de l'affichage de la transmittance. Le programme 
indique la largeur de l'écran exprimé en microsecondes ; il calcule éga- 
lement le temps d'établissement à —-5 %. 

La hauteur de l'échelon étant de 100 pixels, on peut ainsi mesurer direc- 
tement sur l'écran les divers dépassements : un pixel égal 1 % du dépas- 
sement. 
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LARGEUR BE L'ECRAN 280 mis 











Réponse à un échelon : BUTTERWORTH ORDRE 2 


Qi 























Réponse à un échelon : LEGENDRE ORDRE 3 
4 























: TCHEBYSCHEFF ORDRE 4 


Réponse à un échelon 
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Fi 


PONS 














BESSEL ORDRE 5 


. 
. 


Réponse à un échelon 
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10 BEM RRE#RERERÉÉNERERNEFÉERÉRÉÉÉERSE 

2û 

30 

40 

SÛ 

ëû RE 

F6 CLS 

80 SPINT:PRINT" FILTRES ACTITE" 

QÙ PRINT : 

100 FRINT:PRINT" Nous vous rappelons que les principauxtupes de filtr 
es sont!" 

120 PRINT:FRINT"-{ EUTTERWORTH caracterises par une * 

130 PRINT'courbe de reponse plate à l'oriçine." 

140 PRINT:PRINT*-2 LEGENDRE caracterises par une cassure" 
150 FRINT'nette a ia frequence de coupure.”* 

160 PRINT:PRINT®-3 TCHEEVSCEFF caracterises par ur gain” 
170 PRINT'a peu pres constant dans la totalite” 

180 PRINT" de la bande passante," 

190 PRINT:PRINT"-4 BEGGEL optimises cour la requiarite‘ 


+ 


200 PRINT‘ du tençs de pronsgation ce grouse dans" 

PRINT'la bande passsante,* 

7 PRINT:PRINT" TAPEZ V nr CHGIX SVP, puis ENTER" :: INPUT 7 
240 ÏF Z£i OR Zi4 THEN 

250 CLS 

280 PRINT: PRINT'DONNEZ DU FILTRE QUE VOUS" 

2? VOIER { de Z a 5 }.*::INPIJT Ÿ 

290 IF YEZ OR Vis 
300 R=Y:Y=Y-{ 

gi CLS 

320 D=10*G.02:x=1/(100)) 

330 DIM A(200ÿ:DIX B(200: 

gg PLOT 0, JC: DRANR SAR = 3,3: 
330 PLOT 20,10: DRAUR OC, ASC:DREUR -5,-5 
360 FOR N=1 TG 3:FLOT 20+10C4N, 16:DRAU 8 
370 FOR N=0 TO 5:PLOT 20, 160440: ERAUR 8,C:h 
300 ON Z GOSUE 390,865, 1330, 1800 

J90 A$=EUTTERMORTH ORDRE" 

809 ON Ÿ G07G 410,51C,820,740 


1 
2 


hi 
Les 
Les | 





3: DRAUR -5,-5 
3: DRA4R 3,- 


Ê 

d: 

C4 
Ur 
NEXT 
EXT 


DRAWR 
DRAWR 
# 0,8: 


410 REM EUTTERWNORTH ORDRE 2 
32û A 
436 GÜSUE 2 


440 GOSUE 
450 FOR N=! Tü 200 


pas 


480 X=X+D 

470 A(NI=FNS(X):BCN)=FNAIT) 

480 PLÔT 20+24N, 32040041 GG(ACN)) 2,3 
490 NEXT N 

500 6070 2450 

510 REM EUTTERKORTH ORDRE 3 

520 A=i:B=1:C=1:A$=AS+STRE(Y+1) 

520 GOSUR 2340 

540 GOSUR 2270 

550 GOSUR 2290 

560 FÜR N=1 TG 205 

570 X=X+D 

580 A(N)=(FNS(X)#FNPiX))2E(N)= ie ss 
590 PLOT 20+24N,320+80#L0OG{A(N))/2, 
60G NEXT N 

6iG GOTO 2450 

620 REM EUTTERWORTH ORDRE à 

630 A=1:2=1.0477:C=i:A$=AS+CTRE(V+ I; 
640 E=1:F=0,7653 

650 GOSUE 2310 

66G GÜCUR 2290 

670 GOSUR 2330 

680 FÜR N=1 TO 200 

699 X=X+D 

700 A(N)=(FNGEXI#FNOUX)): BÜNIZENACX)+ENEL X) 
710 PLOT 20+24N,320+80#i GG{AIN) )/2,3 
720 NEXT N 

730 GGTO 2450 

740 REM EUTTERUGRTH Sa E 5 

750 A=1:E=1,6i8:C=1:AS=AS+STRS(Y+i) 

760 E=1:F=0.613 

770 GOSUR Z310:GGSUE 2270 

780 GO0SUR 2290 

790 GOSUE 2330 

800 FOR N=1 TG 200 

810 X=X+D 

820 A(NI={FNS(X)#F NQ OO #FNP OC) :BNIA FNAX)+FNECXI4FRCEX) 
830 PLOT 20+24N,320+80#LO0GCA(N))/2,3 

849 NEXT N 

850 GOTG 2450 

BG A$="LEGENDRE GRDRE * 

870 ON Y GüTG 080,960,1090, i£ig 

880 REM LEGENDPE ORDRE 2 
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890 A=i:P=ZA0.5rAË=ASECTREE VE) 
300 cüs VE Re. 
SiG 6 296 

LA FOR N=1 TG 20 


GJC X=Y4L 
940 ACNI=ENGEX REIN 
955 PLGT 
Gé0 NEXT À 
970 GOTO 2250 
Q8G REM LEGENDRE ORDRE 3 
900 A=1,0744:E20, 7417: C2; , 612: AS=AS<STES I Y+i; 
i0GG GÜSUR 2940 
1010 GÉSUR 2790 

C20 GÜSUR 2770 
1030 FGR N=1 TG ZCG 

1040 X=X40 

1050 ACNI=FNSEX 4PNBCX Ds EEN)=ENACX AENCI x: 
1060 PLOT 20+424N, 7204808 0G{AIN):/2,3 

1070 NEXT N 
1980 GGTE 24 
1090 REM LE 
1100 A=2, 32 
1110 E=1.055: 
1120 GOSUE 
11930 GÜSUE 
1140 GGSUR 
1150 FOR N= 
1160 Y=X+0 
1170 AUNI=FNSOX)#FAGX) 2 BNI=ENAUXI HENRI 1} 
1190 PLOT ZO424N, 32D+B04LOGtAN) 72,3 
1190 NEXT N 
1200 GüTG 2450 
1210 REM LEGENDRE GSDSE 5 

1220 A=2,.0115:121,56j4:C=2. 176: At=AR4CTRSIY+1) 
1210 E=1.0406:F29,3196 
1240 GOSUR ZT 
1250 GÜSUR 2290 

F 









es ho 


TG 295 


2 
1260 GÜSUR 2310:G0SU8 7275 

1270 FÜR N=1 T0 290 

1290 Y=X+D 

1290 ACN)=FHSCXIKPNE CE E#PNO NX : BONI=PNACX PNA X HE NCUE: 
306 PLOT 20424, 3204054106 (41 

310 HEXT N 





4 
i 
{ 
i 


1320 T6 2450 

1330 A$="TCHEEYSCHEFF GRDRE * 

1340 ON Ÿ GOT0 2459,1450,:5e0, 1890 
1950 REM TCHERYSCAHECF ORDAE 2 

1360 A=0. 6595 :Ezû, 040: At=AE4CTRS I Y+l 
1370 GOSUR 7330 
1380 GOSUR 2200 
1390 FOR N=1 TG ZCO 

1400 X=x+#0 

1410 A(N)=ENG(X):BENI=ENAUX) 

1420 PLOT ZOHS4N, 120485 CS CAIN) 2,3 
1430 NEXT N 

1440 GÜTG 2:50 

1450 REM TCHERYSCHEFF CEDRE 7 
1460 À=0,8753:5-0,5484:C=i,55é 
1479 GOSUE 2330 

1480 GÔSUR 7290 

1496 GOSUE 7270 

1500 FOR N=1 Tü 200 

1519 X=X4D 

1525 ACNI=ENGOY)#FNCE Ke RIN)=FNAIXI+ENC! 5 
1530 PLOT 20424, 320+804L05{A(N):/2,3 

1540 NEXT N 

1550 GTÈ 2450 

1560 PEM TCHERYSCHEFF GRDRZ à 

1570 A=0, 9402: 125, 7207: A$=h#iGTaéiY+i) 
1589 E=2,8097:F=2,3755 

1590 GSSUR 2310 

1600 SOSUR 2290 

1610 GÜSUR 2330 

1520 FOR N=i TG 200 

1635 X=X+D 

1640 AINI=ENS ON #E NON EUN) NAI Ki+ENRi 
1650 PLOT 20424, 3720460400 (AN) 3/2 
1666 NEXT N 

1670 GCTO 2450 

1680 RE" TCHERYSCHEFF ORDRE 5 

1695 A=2,0974:2=1,2298:t 

1700 E=0.9854:F=0.7151 
1710 GÜSUR 2910 

1720 GOSUE ZZ9C 

1720 GOSUR 2J{0:GOSUR 2275 
1740 FÜR N=! TO 200 
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LAËSAG4ET 28 (Y4L) 











STE 


C=2, 750: AË=AS+STRE VE) 






1750 X=Y#D 





iFëG À LI=ENACKIHFNEE XI HENCE SE) 
1770 PLOT 20+24N, 326 }/2,3 





1780 NEXT N 

1790 GÜTÜ 2450 

1E0Û A$="HECGEL ORDRE” 

1810 ù Y GOT 1820,1920,2030, 3155 
1820 REM DESSEL ORDRE 7 

1820 A=0.618:E=1. 2616: A$=AS4CTRE (Ve) 
1840 GOSUR 2330 

1950 GOSUR 2290 

1860 FOR N=t TO 700 

1870 X=X#D 

1880 AÏNj=FNS(X): Fit 
1895 PLOT DOI 22 ECAL 06 A(Ni)/2,3 

1906 NEXT N 

1910 GOT 2450 

1920 REM REGGEL CADRE 3 

19930 A=0.4771:2=0,9996:C=0,75L:49=A$+CT#siV+i) 
1940 GÜSUR 2330 

1950 GÜSUR 227C 

1960 GOSUR 7290 

1970 FÜR N=1 TO 200 

1989 X=X+D 

1990 AUNIZ(FNS OK EPNS OX) BRENIZFNAUXIHENCEE) 
2000 PLOT ZO+24N, J2D+8C#LOG(AUN))/2,3 

2010 NEXT N 

2020 GGTO 2450 

2040 REM HEGSEL ORDRE 4 

2040 A=G. 3869: R=0, 7742: A$=A#+STRE(Y+ ) 

205G E=0,.49899:F=1,3396 

2080 GÜSUE 3310 

2076 GÜSUR 7290 

2080 GOSUR 2330 

2090 FOR N=i TG 200 

2100 X=X+0 

2110 AÏNIZIFNSEKIAFNQUX 2 BUN)=ENACE)+FNECX) 
2120 PLOT ZOH4Z4N, 320+90#L06{A1N))/2,3 

2198 NEXT N 

2146 GOTG 2450 

2150 REM HEGSSEL CARE 5 

2160 4=0,.4128:E=1, 140 D, ES AS=AEESTREIV+ 1) 
2170 E=9.4245:F=6 mi 


ih 





2180 GOSUR 2310:G0SUR 2270 

2190 GOSUR 2290 

2200 GOSUR 2330 

2210 FÜR N=1 T9 200 

2220 X=X+D 

2230 A(N)=(ENSEX AE NQ OU) #FNPEX) D: BEN =FNACÉ)+FNEE XI HENC(X) 
2240 PLOT 20+24N, 320+804L0G{A(NI)/2,3 

2250 NEXT N 

2260 GÜTG 2450 

2270 DEF FNP(X)=1/({#(C#X)AZ)4G,5 

2200 DEF FNC{X)=ATN(CEX) : RETURN 

2290 DEF FNS(X}=1/{{1-ARXAZ)AG+ (REX) AZ) AQ,S 
29300 DEF FNACX)=ATNE(A#X-1/X3/E) : RETURN 
2310 DEF FNG(X)=1/((1-ERXAZ) ASE (FX AZ) AG, 5 
29320 DEF FNB(X)=ATNE(E#X-1/X)/F3: RETURN 
2330 LOCATE 20,24:PRINT A$ 

29350 LOCATE J,3:PRINT"+10 de" 

2360 LOCATE J3,8:PRINT"-10 de' 

2370 LOCATE J,f1:PRINT*-29 cB* 

2380 LOCATE 3, 1G:PRINT"-30 dh* 

2390 LOCATE J,ié:PRINT*-40 dP* 

2400 LOCATE 6,4:PRINT*W9/10* 

2410 LOCATE 14,4:PRINT'WO* 

2420 LOCATE 19,4:PRINT"w0*x10" 

2440 RETURN 

2450 LOCATE 18,1: PRINT'REPONSE TRANSITOIRE* 
2460 LOCATE 18,2:PRINT" APPUYEZ SUR © + ENTER ss INPUT 29 
2470 CLS:REM REPONSE TRANSITOIRE 

2480 PRINT:PRINT'REPONSE TRANSITGIRE" 

2490 PRINT: INPUT "Frequence de coupure a -3 cE*:ft 
2500 DIM C{200):DIM D{200) 

2510 FOR K=1 T0 30 STEP Z 

2320 C(K)=SO#LOG(ZS#K)/2.3 

2530 NEXT K 

2540 CLS 

2550 T=4/(F0+400) 





2560 PLGT 10,100:DRAUR 470,C:DRANR -5,5:DRAUR 5 
2570 PLOT 10,300:DRAWR 470,0: DRAUR -5,5: DRAUR 5 

2380 PLOT 19,80:DRANR 0,350: DRAUR 3,-5:DRAWR -5 
2600 LOCATE 1,23:PRINT'EARGEUR DE L'ECRAN="; CUT 


“micro-s." 

2620 L=9#PI/40:TR=0:E=0 

2630 IF R=2 OR R=3 THEN P=PI/Z 
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2640 IF R=4 OR R=5 THEN P=PI 
2656 co N=0 TG 199 
2660 U=O:F=Fü/i0 


2670 FOR K=1 TO 30 STE? 2 
2880 U=U(AUCEKD3/KDRBINE (24 TAN ARTE CUKD)-P 






2700 PLOT 10+24N, 200400 U/P: 
2HO IF UC AND E=D THEN TRSTRT 
2720 IF U>L THEN Ezi 

2730 NEXT N 

2740 LOCATE 8,1:PRINT'TENRS DE RE 
2 5 


2750 GÈTO 27 


m 


PONGE :" 


Der al 





Chapitre III 


FOURIER 


SÉRIES DE FOURIER 


Le théorème de Fourier a un grand nombre d’applications..Il a surtout 
été employé comme base théorique de démonstration. S'il n’a guère 
été utilisé en pratique, c’est qu'il conduit à des calculs inextricables. 
Aujourd'hui, l'ordinateur permet de l'exploiter. 

Ce théorème dit qu'à toute fonction périodique, peut se substituer une 
somme de fonctions sinusoïdales, dont les fréquences sont multiples 
de la fréquence fondamentale. Le nombre de termes de cette somme 
est infini ; voici trois exemples de développements : | 





= AA Teese les then Bot 1 
T 3 5 





y= (Sin ut- Sin 3 wt+ Sin DURS" 2 





Figure 3 


_ Sin 2 ut je Sin 3 et 


2A 
Ve 2 3 


LE PROGRAMME N° 17 


Ilest démonstratif. || propose de montrer sur l'écran la forme du signal 
obtenu quand on emploie K harmoniques. On visualise une période et 
demie du signat 3 pris en exemple. Après chaque balayage d'écran, un 
nouveau signal, augmenté d’un harmonique, vient remplacer le précé- 
dent. 


LES VARIABLES 


A(200) : contient la position de tous les points de la courbe sur l'écran. 


K : nombre d’harmoniques pris en compte. 

T : variable temps. La période du signal correspond à 120 points ; 
T varie de O à 179 ; une période et demie du signal est ainsi 
affichée. 

Y : valeur de l'harmonique K au temps T ; A (T) est la somme, 


au temps T, de tous les harmoniques de rang inférieur à K. 
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RECONSTITUTION D'UN SIGNAL 
EN DENTS DE SCIE AVEC : 


1 terme (le fondamental) 


doulez-vous superposer les courbes?70/H?7 


K= 1 


Da _ 


2 termes (le fondamental et le premier harmonique). 


1 REM RIRE 

20 REM # # 

30 REM # SERIE DE FOURIER * 

40 RE * * 

50 REM HEIN 

ë0 ? 

70 INK G,O:INK 1,6: PAPER G:PEN 1 
8G DIM 4(200) 


90 CLS:PRINT "Voulez-vous suerco s (EN)? 


n 
mn 
ss 
+ 

vm 
in 
nm 
[a] 
a 

+ 
ue 

10 
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SERIES DE FOURNIER 


DONFeE VON superposer les courbes?70/N? 


ESS K= 2 














S'uHes voue superposer les courbes?70/N? 





K= 5 


ses 





DONTER Ve superposer les courbes?7O0/N? 


A] 
Us. 


K= 15 





95 A$=INKEYS:IF 4$=*" THEN 95 
CG IF At="0" THEN N=i 

iii CES 

20 FOR K= 1 TG JO:REM K EST LE NOMBRE D'HARMONIQUES PRIS EN COMPTE 
150 L=1 

10 IF INT(K/2)=K/2 THEN L=-1 
186 E$="K="4STR#iK) 

210 LOCATE Z8,5:PRINT B$ 

240 FOR T=C T0 175 

250 Y=L#SIN(K#PI#T/60:/K 

270 A(T)=A(T)+Y 

310 PLOT 347, 200-1004 (T) 

J20 NET T 

330 :F K=30 THEN END 

335 IF N=1 THEN 430 

340 €LS 

350 NEXT K 


APPLICATION DE LA SERIE DE FOURIER 


Lorsqu'un signal non sinusoïdal est appliqué à un circuit électrique, il, 

peut être difficile de déterminer le régime permanent correspondant à 

la tension de sortie de ce circuit. Nous allons voir qu’un ordinateur, dont 

la capacité de calcul est infinie, permet de résoudre un tel problème. 

Nous proposons donc : 

— de relever (ou de calculer) le module et la phase de la transmittance 
du circuit étudié en régime sinusoïdal permanent, ce qui en général 
est aisé à faire ; 

— de construire le signal de sortie du circuit, en le remplaçant par sa 
série de Fourier, chaque terme de celle-ci étant ‘‘pondéré‘’ par le 
module et la phase de la transmittance du circuit à la fréquence cor- 
respondante. 


LE PROGRAMME est très simple. Il prend comme exemple le cas d’un 
circuit intégrateur RC du premier ordre, soumis à une tension 
rectangulaire. 





FRERES @ = — Arc tg RCw 


RC est la constante du temps du circuit, w la pulsation du signal sinu- 
soïdal incident. 

Ce programme permet de comparer sur l'écran la réponse calculée par 
la série de Fourier du signal d'entrée et la réponse exacte, définie par 
une équation du type : 1 -exp{— T/RC) (en toute rigueur, cette réponse 
n'est ‘’exacte’’ que dans la mesure où RC est le dixième de la période 
du signal). 

Remarquons toutefois que le cas pris en exemple est un peu particu- 
lier, dans la mesure où il s’agit d’un filtre de type passe bas et donc 
pour lequel les harmoniques de rang élevé ne jouent aucun rôle. 

La période du signal To est égale à 200 (points). 


Quel est l'intérêt de cette méthode ? En effet, on pourra constater l’ex- 
trême lenteur de celle-ci en faisant tourner le programme. Le principal 
avantage que l’on peut en tirer, et qui est loin d’être négligeable, est 
que cet algorithme fournit directement le régime ‘‘permanent”’ du cir- 
cuit, sans passer par le régime transitoire. Le tracé direct de la solution 
(en bas sur les photos) ne conduit à la solution que si RC«To. Pour 
RC>To, il faudrait tenir compte du fait que la valeur moyenne de la 
tension de sortie est nulle (signal aligné sur l'axe). Même une intégra- 
tion directe de l'équation différentielle fournirait, dans un premier temps, 
le régime transitoire du circuit. 

Nous utiliserons cette méthode dans le programme consacré aux fil- 
tres du 2e au 5° ordre. 


LE PROGRAMME 


Il suffit d'introduire, sur la demande de l'ordinateur, la valeur du terme 
Q=RC/To, seul paramètre jouant sur la forme des courbes. 


Ligne 90 : définition de la solution ‘‘exacte. 
Lignes 130 à 160 : calcul de la tension de sortie pour la série de Fou- 
rier. 


Lignes 170 à 200 : affichage des points et des axes. 


LES VARIABLES 


A et B sont les paramètres de calcul RC/To=0Q. 
A=27RC/To ; B=2#x/To ; ici, To = 200. 


K : numéro de l'harmonique calculé. 
YK _: harmonique de rang K. 
Y __: fonction calculée avec K harmoniques. 


10 REM HAMMEE EE 

20 REM * 4 : = 
30 REM # UTILISATION DE FOURIER +# É ed 

40 REM * ; he 

SC REM HOMME ENEEOEUEE à F 

60 CLS TT —-— 

70 INPUT "DONNEZ LE RAPPORT RC/T DESIRE" :û 

BO A=2*PI%0:P=PI/100 


85 CLS:PRINT *RC/T=":Q b PSE 
90 DEF FNZ(T)=PIx({-EXP(-T/(20040)))/4 © 4 

100 INK O,D:INK 1,26:PAPER O:PEN 1 5. ; 

110 FOR T=0 TO 200 ns 

120 Y=0 É 


130 FOR K=1 TO 31 STEP 3 

140 YK=C1/ (KE (1 (KAA) AZI AC, SISRGINIKERET-ATNEKEA) ) 

150 Y=Y+VK 

160 NEXT K 

170 IF Té=99 THEN Z=FNZUT):PLOT 247+20, 10041207 

180 IF T:09 THEN Z=FNZIT-100):Vi=(éOxPI/4)-60#7:PLOÔT Z4T+2C, (00+24vi 


190 PLÔT 247420, 100+t804P1/4) 
200 PLOT 2*7420,300+604V:PLOT 247420, 200 
210 NEXT T | 
220 END 
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UTILISATION DE LA SERIE DE FOURIER 


Le programme précédent a montré la méthode qui permet de trouver 
le régime permanent d'un circuit dans le cas d'une sollicitation quel- 
conque, dans la mesure où l’on connaît à la fois le développement en 
série de Fourier de la sollicitation et la réponse en fréquence (module 
et phase) du circuit considéré. Nous allons utiliser cette méthode pour 
trouver la réponse d'un filtre passe bas du second ordre à un échelon 
(problème traité par ailleurs dans cet ouvrage, puisque dans ce cas par- 
ticulier, on connaît l'équation différentielle qui régit le fonctionnement 
d'un tel circuit). 

Le circuit étudié est donc le suivant : 








— Vs 1 
FO 1+2m +6 
wO wO 
1 _ Lwo 


; 1 
ec Q = —— _ 
_N mn RO TT 


c'est-à-dire, en module et en phase : 





IT ol = =" 
V1 2 + (2m 272 
&wO? &wO 
; @w= Arg (T)= — Arc tg “em 9/00 
1 —-w?/w?20 


Chaque tenue du développement en série de Fourier de Ve, de la forme 
1/K sin K w t devient donc, en sortie : 1/K x ITul x sin (Kwt — pu). 


La recomposition de la série de Fourier après traitement de chacun de 
ses termes, donne la série de Fourier relative à la tension de sortie. 
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LE PROGRAMME 


On pose SE =C, en appelant wF la pulsation du terme fondamen- 


wO 
tal. Ainsi, pour chaque harmonique de rang K : 


© 2 Ko _ KC. 
wO wO 


Nous étudions ainsi un signal rectangulaire, dont la fréquence est C fois 
la fréquence propre du circuit. D'autre part, pour le fondamental (K= 1), 
on a: 


@Ft = 2at 
200 





ce qui fait que 200 points sur l’écran correspondent à une période, aussi 
un choix est-il possible au début du programme : affichage d'une demi- 
période (réponse à un échelon) ou d'une période complète (réponse à 
un signal rectangulaire). Comme pour le programme précédent, nous 
obtenons dans tous les cas le régime permanent (le régime transitoire 
est étudié par ailleurs dans ce livre). 


Lignes 130 à 210 : calcul de la tension de sortie par le développement 
en série de Fourier. 


LES VARIABLES 





m  : coefficient d'amortissement du circuit. 
T : variable temps. 
K _: rang de l'harmonique considéré. 
2 
A : valeur de 1 -£ 
w?20 
B 


: valeur de 2m w/w0 pour le rang K. 

AK  : inverse de l'amplitude de l'harmonique de rang K. 
PH _: valeur de okw. 

YK _: valeur de l'harmonique de rang K au temps t. 

Y : somme des harmoniques. 


O ROM RHAHRAIRRARRERRERRREREERERÉÉRÉE 

C REM +# L 

Q REM # UTILISATION DE FOURIER {Ji +# 

40 REM * 4 

SC REM HIER EE PRET 

&0 TO=200 

70 CLS 

80 INPUT "DONNEZ LE COEFFICIENT D'AMGRTISSEMENT 
CL] :M 

GG PRINT: INPUT "DONNEZ LE RAPPORT Fi/Fo*sC 


Ci RO va 


_ 
DE 


100 PRINT:INPUT VOULEZ-VOUS LA REPONSE À UN ECHELON 


RECTANGULAIRE (2) " ;N 
110 CLS 

120 1=0 

130 FOR T=0 TO 10CAN STEP N/2 

140 FOR K=t TO 10/C STES 2 

150 A=1-{KkC)°2 

160 P=ZaMkEC 

170 AKSHRL{A"2+222)40,5) 

180 PH=ATN(R/4) 

190 IF PHIO THEN PH=PI+Ph 

200 VKECL/AK)RSINC(Z4K#PI4T/200) -PH) 
210 Y=v+v 

220 NEXT K 

220 PLOT 20+2#1,200+100#Y 

240 PLOT 204281, 200 

250 Y=0 

260 1=1+1 

270 NEXT T 

280 END 


VOTRE CIRCUIT 


{150U À UN SIGNAL 


Remarquer que le choix du nombre d’harmoniques pris en compte est 
asservi au rapport C =Fi/Fo. En effet, on peut estimer que tous les har- 
moniques de fréquence supérieure à 10 Fo sont complètement élimi- 
nés par le filtre (atténuation de l’ordre de — 40 dB, soit atténuation d'un 
rapport 100). Le dernier harmonique à considérer a donc une fréquence 
de l'ordre de 10 Fo, d'où, pour la plus grande valeur de K : 


KFi= 10 Fo=KCFO—KC = 10-K= © 
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CALCUL DES COEFFICIENTS 
D'UNE SERIE DE FOURIER 


La série de Fourier d’une fonction périodique s'écrit : 


2rt 


TT + A2 sin (2x ER 








f(t) = Ao + A1 sin 


2xt 27t 


+B2 cos (2x ) +... 








+ B1 cos 


Les coefficients se calculent par les relations : 


T io dus An 2 À ftixsininx 22 jdt 
0 7e T 


et Bn= oh f(t) x cos{n x 2rt jt 
T . T 


Le calcul n’est aisé que pour des fonctions f{t) simples. 
Nous proposons de remplacer les intégrales précédentes par des som- 
mes en nombre fini. Ainsi : 





A0 


I 
-|- 





— f{t) x sin(n x £= )dt devient 


2xt 
T 


t variant de O à T par bonds de At:At=1 (point). 

Nous définissons la période de f(t) comme étant égale à 200 (points). 
Un tableau contient les 200 valeurs prises par f(t) sur cet intervalle. 
Nous faisons le calcul dans le cas d’un signal rectangulaire (f(t) = 1 pour 
1<t<100 et —1 pour 101<t=< 200). La série de Fourier correspon- 
dante est alors connue et permet de chiffre la précision de la méthode 
(6% sur 21° harmonique). Les termes trop petits et non significatifs 
{inférieurs à 10:4) sont éliminés. Une meilleure précision pourrait être 
obtenue en définissant f{t) sur un plus grand nombre de points. 


27 : 
— Z f{t) xsin(n x 
T o 





)x At 


LE PROGRAMME 


Il y a un affichage par ligne sur l'écran du numéro de l'harmonique et 
des valeurs de An et de Bn. 


Lignes 120 à 190 : calcul de la somme correspondant à l'intégrale à 
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évaluer pour l'harmonique de rang N. 

Ligne 210 : calcul de la valeur exacte de l'amplitude de l’har- 
monique de rang N (seulement valable pour 
l'exemple proposé, à titre de comparaison. Sup- 
primer cette ligne pour d’autres cas). 

Ligne 240 : mise en mémoire des valeurs calculées. 

Lignes 270 à 290 : définition de la fonction étudiée. Peut-être rempla- 
cée par une autre fonction, de période égale à 200 
points. 

Lignes 310 à 340 : mise en page sur l'écran. 

Lignes 350 à 390 : mise en page sur l'imprimante. 


LES VARIABLES 
A(200) : valeurs de la fonction étudiée [f{t)]. 


S(30) et C(30) : matrices contenant les valeurs calculées des coeffi- 
cients An et Bn. 


E(30) : contient les coefficients (An) exacts de la série des 
sinus. 

N : numéro de l’harmonique. 

T : variable de temps. 

ST : valeur d'un élément de la somme correspondant au 
calcul de An. 

CT : idem, mais pour B. 

SE : valeur exacte du coefficient An. 

U : variable de boucle locale. 

N SINUS COZSINUS 

calcule exact calcule exact 

i 1.275 1.275 — 020 (a 

OO GG GC e 

3 G.4%4 Q.474 SO La 

&  —,000 CO, GO — , DC CG 

nm 0,554 GET — AQE0 Q 

de [SNS QC QC [el 

7 0.iSsi 9.18% —,020 te 
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10 RM RRERRERERERERRREEE 





20 REM # + 
30 REM # 4 
4Q REM + D° £ 
50 REN *# 4 
60 REM # £ 
70 REM HRRRRRARRAEENREEE 





80 CLS:INPUT "DÉSIREZ VOUS UNE 
100 DIM A{200):DIM S(30j:DIM C{ 
110 GOSUR 270 


4 


120 FOR N=i TÜ 30 





136 S=C:C=0 

14G FOR Tzi Ti 200 

150 ST=AÎT)#CINUZ#PIKN#T/200) 
1éG S=C+HÛT 


170 CT=ATI#COSC2APTRNET/200) 
180 C=C+CT 

190 NEXT IT 

200 S=6/100:C=C/160 


210 IF INTIN/2)-N/220 THEN SE=4/ (PIAN) :GE=INT(SE#10000/19000 £ 





= 
220 IF A$="0" THEN GOSUR 380 

230 PRINT Ni: PRINT TAB(S);:PRINT USING". H#8° 36; : PRINT 
USING"A HA": GE; :PRINT TAB(21)::PRINT USING'#.###" 305: 


Z4G SENI=S: CINI=C:EN)=SE 

250 NEXT N 

260 END 

270 REM EXEMPLE SIGNAL CARRE 

280 FOR U=0 Tü 10G:A(U)=t:NEXT 
290 FOR U=101 Tü Z0C:A(U)=-1:NEXT 
300 CLS 

310 PRIAT® N° TAR(9) "SINUC* TAR(23) "COSINUS" 

320 PRINT TAR(S) “calcule ekact calcule evact":PRINT 
330 IF A$="0" THEN GOSUE 350 

340 RETURN 

430 PRINT #6," N° TARI9) "SINUS* TAB(2 
JÉ0 PRINT gg," TARES) ‘calcule exact calcule exact* 
D RETURN 


AUjC® 


NV 





SPRINT #0, USING®H HAN IGES: PRINT #8, TAR(Z1) sp 
HER ICI RINT #9, TAB(Z9) °° 
390 RETURN 


INT ee 


D PRINT RON TARES) "ss FRINT #9, USINGH.###" SPRINT # 


PRINT 
B{29}:"C 
TAZ{13 
USINE E "4 


Chapitre IV 


RÉGIMES TRANSITOIRES 


REGIMES TRANSITOIRES 
DU PREMIER ORDRE 


Le programme propose une étude systématique des circuits RC du pre- 
mier ordre, pouvant être soumis à différentes sollicitations (signaux 
triangulaires, rectangulaires et sinusoïdaux). Le régime étudié est donc 
le régime transitoire de ces circuits. En laissant le programme se dérouler 
un certain temps, on voit apparaître le régime permanent. La gamme 
des signaux proposés n’est pas limitative. On peut en utiliser d'autres 
qui seront définis, soit sous forme de fonctions, s'ils ont une expres- 
sion analytique FNA(T), soit sous forme de d'éléments de tableau A(T). 
Le problème est résolu par l'intégration numérique de l'équation diffé- 
rentielle. L'étude de tout circuit défini par le même genre d’équation 
se fera donc de manière identique. 

Le programme se présente sous forme de menu. Il faut d’abord choisir 
le circuit dont on veut étudier le fonctionnement (tableau 1), puis la 
forme du signal et leurs caractéristiques (tableau 2). Ensuite, l'écran 
montre le signal d'entrée choisi et le signal de sortie correspondant. 

Comme toujours, les premières courbes obtenues sur l'écran représen- 
tent le régime transitoire du circuit. En général, et sauf cas particulier, 
le second passage représente déjà le régime permanent, c'est-à-dire 
celui que l’on observerait sur l’écran d'un oscilloscope dans les mêmes 
conditions. 


LE PROGRAMME 

Lignes 90 à 190 : définition des dessins 

Lignes 900 à 1540  : premier tableau : choix du circuit 
Ligne 1560 : deuxième tableau : choix du signal 
Ligne 1610 : introduction des valeurs numériques 


Lignes 1700 à 1750 : initialisation des calculs 
Lignes 1760 à 1910 : caculs : 


Lignes 1990 à 2020 : sous-programmes relatifs à chaque circuit 
Lignes 2030 à 2110 : rappel des paramètres en cours sur l'écran 


LES VARIABLES 


N : numéro du circuit choisi. 

K : numéro du signal choisi. 

R, R1, R2, C : valeurs numériques des éléments du circuit choisi. 
F : fréquence du signal d'entrée. 


vo : tension initiale (à t=0o) sur le condensateur. 
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: tension d'entrée (en V) crête à crête (elle est supposée 


sans composante continue). 


la discontinuité des signaux 


DT : intervalle de calcul, pris égal à 1/(200 F). L'écran 
complet contient donc une période du signal d'entrée. 

T : variable temps (il ÿ a 200 intervalles sur une période). 

L : variable imposée par 
considérés : L=0 pour la première demi-période, L=1 
pour la seconde demi-période. 

A, D : variables de position relatives à VE (tension d'entrée : 
valeur crête) et à VS (tension de sortie) pour l'impri- 
mante). 

VE : valeur instantanée de la tension d'entrée. 

| : intensité instantanée du courant dans le circuit. 

DV : Variation de la tension aux bornes du condensateur pen- 
dant l'intervalle de temps DT. 

VC : tension aux bornes du condensateur à l'instant T. 

FNA(T) : fonction correspondant à la tension d'entrée. 


Pour la partie du programme comprise e 
rera à l'annexe relative à l'imprimante. 


ntre 2000 et 3050, on se réfé- 








F Cac Ki © KL 
+ CR nRE Re 
4 Li ; , 
FF Li 1} LL 
Ds re Te 
£ res a 4. 
CHOISISSEZ LB CIRCEIET GOMT LOIR 
LIEN RMEULIER LE FONCTIONNEMENT 
FE 
Œuel signal voulez-vous lui appliquer? 
Æ Kegtanaquiaire 13 
HF ÉXIRMNMSGUIRIrE {223 
He Sinaissoidalie LE 
ZE 
Guelles sont es caracteristiques 
HZ LÉSESSIE 
C=+ 1E-S 
Freuence dus signal d'entresT 15645 
Tension d'entreer 18 
Tension initiales 2=ux ia capwaciter 6 
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(1) Signal rectangulaire 
V=10 V crête à crête 
Vo=0 V 


ACL FE 
ee 1666 HE 





(2) Identique à 1 
mais Vo= —-10 V 


CIRCUIT Ho 1 

E= 16668 

C= S.6H0BB001L F 
Frequencez 16 HE 





(3V=10Vc.àc. 
Vo=0 V 





(4) signal rectangulaire 
V=10 V 
Vo=0 V 


CIRCUIT Mae 2 
HE 16698 
L= 5.698668 81: 


Frequence— So hz 
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(5) V=10 V 
(suite du dessin 4) .Vo=0 V 





(6) V=10 V 
Vo=0 V 


LE UT bel 
% 


fe dE 
rt Le 


tu Li 


TE 
% 4 


th 
4 


ca pra 
DRE IE TE| 
th 
de Kid 


fé 
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(8) Signal triangulaire 
V=10 V 
Vo=0 V 





(9) Signal sinusoïdal 
V=10 V 
Vo=0 V 


CLIRCUET Ho à 
E= 1865 

L= 8. ee ie FE 
Frequence— [SISTS 
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(10) Signal sinusoïdal 
V=10 V 
Vo=0 V 


SIS] 
G.G65BE6601 F 
cquencez 1666 Hz 





10 REM HER 
20 REM * 4 

30 REM # REGIMES TRANSITOIRES # 

40 REM# DU PREMIER ORDRE  +# 

50 REM * + 

60 REM HHHHHRENNERENEEAEREREEES 

80 CLS 

90 REM DEFINITIONS GRAPHIQUES 

100 c=150 

105 SYMROL AFTER c 

110 SYMBOL c,0,63, 32, 224, 32,63,0,0 

120 SYMEOL c+1,0,252,4,7,4,252,0,0 

130 SYMPOL c+2,40, 40,40, 239, 40,40, 0 

140 SYMBOL c+3,8,8,62,34,34, 34, 34,34 

150 SYMBOL c+4,34,34,34, 34,34, 62,B,8 

160 SYMBOL c+5,8,8,255,0,255,8,8,8 

170 SYMEOL c+6,0,0,0,255,8,8,8,8 

180 SYMPOL c+7,8,8,8,8,255, 74, 148,0 

190 SYMBOL c+8,0,126,129, 129, 68,26, 231,0 
900 REM DESSIN DES CIRCUITS 

910 LOCATE 1,5:PRINT CHR$(c):CHR#(c+1):CHRS(c+é) 
920 LOCATE 10,5:PRINT CHRS(c+2):CHR$(c+6) 
930 LOCATE 20,5:PRINT CHR$(c) ;CHRS(c+1):CHRS(c+2) :CHR$(c+b) 


940 LOCATE 30,5:PRINT CHR$(c);CHR$(c+1):CHR$(c+é) 
un 


950 LOCATE J,6:PRINT CHR$(c+5) 
960 LGCATE 11,6:PRINT CHR$(c+3) 
970 LOCATE 23,6:PRINT CHR$(c+1: 
980 LOCATE 32,6:PRINT CHR$(c+3) 
990 LOCATE 3,7:PRINT CHR$(c+7} 
1000 LOCATE 11,7:PRINT CHA$(c+4) 
1010 LOCATE 23,7:PRINT CHR$(c+4) 
1020 LOCATE 32,7:PRINT CHR$(c+4) 
1030 LOCATE 1i,8:PRINT CHR$ic+7; 
1040 LOCATE 23,B8:PRINT CHR$(c+7) 
1050 LOCATE 42,8:PRINT CHR${c+5) 
1060 LOCATE 32,9:PRINT CHR$ic+7) 
1300 REM IMPRESSIONS SUR L'ECRAN 
1310 LOCATE 1,4:PRINT °R° 
1320 LOCATE 10,4:PRINT °C* 
1330 LOCATE 19,4:PRINT "RI 
1340 LOCATE J0,4:PRINT “hi 
1350 LOCATE 4,6:PRINT "C * 
1360 LOCATE 12,6:PRINT "R° 
1370 LOCATE 24,6:PRINT “#Z" 
1380 LOCATE SJ, 6:PRINT "R£* 
1390 LOCATE JJ,B:PRINT °C 
1400 LOCATE 2,10:PRINT "1" TAB(IQ)2 TAB(22) 3 TAB (31) 4 
1410 PRINT:PRINT:PRINT CHOISISSEZ LE CIRCUIT DGNT VOUS* 
1420 PRINT:PRINT VOULEZ ETUDIER LE FONCTIONNEMENT" 

1540 INPUT ""iN 

1550 CLS 

1560 PRINT:PRINT "Quel signal voulez-vous lui apoliquer?" 
1570 PRINT " + Rectangulaire (1)° 

1580 PRINT “ + Triangulaire (2)° 

1590 PRINT “ +* Sinusoidale (3° 

1600 INPUT *":K 

1610 PRINT:PRINT:PRINT "Quelles sont ses caracteristiques" 
1620 PRINT:IF N=1 OR N=2 THEN INPUT “R=":R:GOTO 1850 

1630 INPUT "RIZ": R1: INPUT "R£=?:R2 

1640 R=RI+R2 

1650 PRINT:INPUT "C=":C 

1660 PRINT:INPUT "Frequence du signal d’entree";F 

1670 PRINT: INPUT "Tension d'entree'iV 

1680 PRINT:INPUT “Tension initiale sur la capacite*;V0 
1690 CES:DT=1/ (2004F) 

1700 VC=VC 

1710 GOSUR 2O4O0:REM REVENIR ICI POUR UN AUTRE PASSAGE 


Ca] 


1720 L=0 

1730 IF K=1 THEN DEF FNAÎT)=i 

1740 IF K=3 THEN DEF FNAÎTI=SIN(Z#PI#1/100) 

1750 IF K=2 THEN DEF FNAÎT)=-1+7/35 

1760 FOR T=0 TO 45 

1770 PLOT 10+2+7+200%L, 300 

1780 GOSUR 1920 

1790 ON N GOSUE 1950,2000, 2010, 2020 

1800 NEXT 

1810 IF K=1 THEN DEF FNAÎTI=-1:PL0T 110+200#,300:DRAUR D, -160 
1826 IF K=2 THEN DEF FNAÎT)=3-T/55 

1830 FOR T=50 TO 160 

1840 PLOT 19+2+7+20C#, 300 

1850 GOSUR 1929 

1860 ON N GOSUE 1990,2000,2910, 2026 

1870 NEXT 7 

1895 IF K=1 THEN PLOT 210+200%L, 380: DRAUR G,-160 

1890 L=L+1 

1900 IF L=1 THEN GÔTO 1730 

1910 IF L=2 THEN Ci5:G0T0 1740 

1920 VE=U#FNAÎT)/7 

1930 D=FNAÎT)#50 

1940 PLÜT 109 +24T+2004L, JOD+O04FNAÎT) 

1950 DV=DT#(VE-VC)/{R4C) 

1960 VC=VC+DV 

1970 I=iVE-VC)/R 

1980 RETURN 

1990 PLOT 10+24T+200L, 500+1804VC/ 4: A=1G0RVC/V: RETURN 
2000 PLOT 10+247+200%L, 300+160kR& 1 /Vr AZ DCHRET/V RETURN 
2010 PLOT 10+247T+200%1,300+160kR241I/V:A=LOGRREET /V: RETURN 
2020 PLOT 10+2+#7+200#L, 400+160#(R2# 1400) 42: AZ1004 (RTE T4VC) /Ve RETURN 
2040 LOCATE Z,20:PRINT "CIRCUIT No'sN 

2050 IF N=1 OR N=2 THEN PRINT "R=";R;CHRS(c+0):G0T0 2089 
2060 PRINT "RI=":RI;CHR$(c+8) 

2070 PRINT "R2=";:R2:CHR#(c+8) 

2080 PRINT “C=";C;"F" 

2090 PRINT "Frequence=":;F:"Hz* 

2100 RETURN 

2200 RETURN 


RÉGIMES TRANSITOIRES DU SECOND ORDRE 


Ces régimes sont importants en électronique. Ils ont lieu dans tous les 
circuits du type filtre passe-bas dont le schéma de principe est le sui- 
vant : 





1 


OR AT 
VE 1+2jm w/w0 + (jw/w0)? 


avec ici 2 m=RCwo 


m et wO sont des constantes ‘’universelles”’ dans la mesure où tous 
les circuits ayant mêmes valeurs de m et de w0 présenteront, pour une 
sollicitation donnée, la même forme de réponse. L'étude de ces régi- 
mes transitoires, si elle connaît une solution analytique, est assez déli- 
cate. On se contente en général, de tracer quelques points donnés par 
des courbes universelles. 


L'équation différentielle qui régit le fonctionnement de tels circuits est 
la suivante : 


2m , dVs 1, d'vs 


VE=VS + 
wO dt wO? dt? 











dans laquelle m est le coefficient d'amortissement et wo la pulsation 
propre. Le problème est résolu de la façon suivante : 


(1) à l'instant t, introduction de VE et des conditions initiales VS, dVs/dt. 
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(2) on en déduit : 





























2 
d?2Vs = HO (VE - VS 2m . dVs j 
dt? &O dt 
2 
(3) soit d (dVs) _ d'vs x dt 
dt dt? 

(4) puis dvs = (ancienne valeur de NS )+d{ dvs }) 
dt dt 

(5) enfin VS = (ancienne valeur de VS)+({ ae x dt) 


(6) Ces nouvelles valeurs de Vs et de sont alors introduites 





dans (2) pour permettre une nouvelle étape de calcul. 


Remarquons qu'au départ, il faut introduire des conditions initiales (Vs 
et LS à l'instant t=0) pour permettre le démarrage du calcul. 





La précision des résultats obtenus est liée au choix de l'intervalle de 
calcul dt. Il faut que dt soit ‘’petit’’ devant 27/w0 (période propre du 
circuit), et que d'autre part Ve varie assez peu pendant le temps dt. 
Nous proposons l'étude de ces circuits pour des tensions d'entrée du 
type : échelon, impulsion. Celles-ci pourront être remplacées par d’au- 
tres fonctions. 


LE PROGRAMME 


Il demande d'abord l'introduction des valeurs numériques nécessaires 
au calcul : M et FO. L'intervalle de temps DT est pris égal à 





= , Sauf si l'utilisateur en propose un plus petit. La courbe se des- 
sine ensuite. 

Lignes 90 à 220 : introduction des données. 

Lignes 230 à 260 : initialisation des calculs. 

Lignes 280 à 350 : calculs. 


Sous-programme 380 : graduations, axes, inscriptions sur l'écran. 


LES VARIABLES 


VO : valeur de la tension aux bornes de C. 


DO : valeur de ae ) (au temps t). 





VS : Variation de tension aux bornes du condensateur à l'instant 
t 





D1 : variation de | as ) à l'instant t. 

VE : amplitude de la tension d'entrée. 

M : coefficient d'amortissement du système étudié. 
FO : fréquence propre du circuit étudié. 

W2 : pulsation propre (au carré). 

Wo : pulsation propre. 


A(200) : matrice contenant la fonction d'entrée. 
B(200) : matrice contenant la fonction de sortie. 


N : variable d'incrémentation, pour la définition de la fonction 
d'entrée. 
K : constante de valeur 2M/WO. 
DT : intervalle de calul, pris ici égal au 50e de la période propre 
du circuit. 
2 
D2 : valeur de CR ) au temps t. 





dt? 


10 REM HA  ROOE 
20 REM + # 
30 REM + REPONSE D'UN CIRCUIT # 
40 REM * DU SECOND ORDRE + 


56 REM *# # 

60 REM HIER EEE 

70 CLS 

80 PRINT:PRINT 

90 PRINT" REPONSE D'UN CIRCUITS :PRINT" DU SECOND Q 
RDRE" 


100 PRINT:PRINT:PRINT'FORME DE LA TENSION D’'ENTREE:" 

105 PRINT:PRINT'Echelon (1) eu Impuision (2)":;:INPUT C:IF Céi O8 C2 T 
HEN 70 

119 DIM A(20G):DIM B(200) 

120 IF C=2 THEN FOR N=1 TO Z00:A(N)=0:NEXT N:A(0)=10 

130 PRINT: INPUT "FREQUENCE PROPRE DU CIRCUIT" ;F0 

140 DT=1/(50#F0) 

150 1F C=1 THEN FOR N=0 TO 200 :A(N)=1:NEXT 

160 PRINT: INPUT "COEFFICIENT D'AMORTISGEMENT":M 

170 PRINT:PRINT" Je vous propose d’etudier la courbe dereconse sur: 
JUSINGY#, HHAAAA® : {4/FO}:: PRINT" sec." 

180 PRINT:INPUT" Avez-vous une autre preference(0/N)":A$ 

190 IF A$="0" OR A$="o" THEN PRINT: INPUT * Quel est votre choixlen mi 
cro-s. SVPI":TT 

200 TT=17*0. 000061 

210 IF 17/2060 ÿ DT THEN PRINT'CE N’EST PAS RAISONNARLE....!":PRINT JE 


PREFERE MON CHOIX." 
220 IF (TT/200) €DT AND TT <ÿQ THEN DT=T7/209 
230 VO=0:D0=0:VS=V0:D1=D0 
240 WO=2*PI4FG:H2=W0A2 
250 K=2#M/WD 
260 CLS 
270 GOSUR 380 
280 FOR U=9 TO 200 
290 VE=A(U) 
300 D2=W2#(VE-VO-K#D0) 
310 D1=D2#DT:D0=D9+D! 
20 VS=DO#DT : VO=VO+VS 
330 B(U)=VO:PLÔT 20+2+U, 200+100*B{U) 
340 PLOT 20+2xU, Z00+100VE 
350 NEXT 
369 GOTO 360 
380 REM GRADUATIONS 
390 PLOT 20,200:DRAWR 420,0 
400 FOR N= O TO 200 STEP 50 
410 PLOT 20 + ZxN,196:DRAMR 0,8 
420 NEXT 
430 PLOT 440,200:DRAWR -5,5:DRAUR 5,-5:DRAUR -5,-5 
440 PLOT 20,80:DRAHR 0, 300:DRAMR 5,-5:DRANR -5,5:DRAUR -5,-5 
460 LOCATE 1,22:PRINT'Ccef. d’amortissement:";M 
470 LOCATE 23,15:PRINT USING'#.H#AAAAE: (SOHDT js: PRINTsec/civ." 
490 RETURN 





REPGSGNSE D'UN CIRCUIT 
Di SECOND ORDRE 
FORME DE LAS TENCION D’'ENMTEHEZ: 
Echelon 13 où ZInpulsion 227 3 
FRERQUENCE FROPRE Di CIACUILTTF i£< 
COEFFICIENT D’'OMONITISCEMENTS 2 
Je vous vcropose d'ebudier la courbe de 

CELONSE Six B.2SE-SS sec. 


831 est votre choix(es misrg-s. SMo)b-z 
à 





Coef, d'amaert 








| B.i9E-Q2sec/4iv,. 


sement: 2 


CT 
((] 











=: 
ï 
i 
: = 
ï # 
: : 
‘ 
; B.1R0E-35Tsec/niv. 


(a 


sef. d'amortissement: 1 


séanete 





G.16E-GZsec/aiv. 


LCoef. d'amortissement: 9.787 








st 














REPSNSE L’UN CIRCUIT 
Di GECOGRD VRDRE 
DT LS TENSION D'ERTERLE: 

Qu Impuilsios LEXT 1 
CGPRE DEL CIRCUITS ES 
D'AMORTISCENENTT .5 
ovose A'etiuiier l4Aa courbes Ge 
S.4GE-H5S sec. 

Avez-vous une autre preferenceco/"H3>z 











S.1BSE-GZTsec/div. 


Coef. t'anortissement: B.5 





menti 





D.1HDE-QZsec/dis. 


Coef. d'amortissement: 1 
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ns ee 
. 


. À  BMNÉE-G3sec/div. 


2 


is 
| 


Coef. d'amortissement: G.8T 


BDB.1SE-9ZT5ec/ûit. 


LCoef. d'amortissement: 9.1 


el 
[+] 
î 
M 


DR ne nn D be 


+ 


B.10E-9Z5sec/diuv. 


R'amortissement: 9.95 





Chapitre V 


CIRCUITS NON LINÉAIRES 


DIODE EN COMMUTATION 


La diode semi-conductrice est un composant essentiel en électronique. 
En commutation, elle présente des caractéristiques qui dépendent de 
sa technologie. Le choix d’un modèle particulier ne peut alors être quel- 
conque. 

Nous allons étudier les formes du courant et de la tension aux bornes 
d'une diode lorsque celle-ci est soumise à une tension rectangulaire. 





Nous allons introduire pour cela un ‘’modèle”’ par parties de la diode, 
faisant intervenir les phénomènes de charge stockée et de capacité de 
transition. 

Nous distinguerons ainsi plusieurs types de fonctionnement : 
1>0:V=0,6 V : dO/dt +Q/T=I|, en appelant Q la charge stockée dans 
la diode et T la durée de vie des porteurs minoritaires de celle-ci. 
1<O et Q>0 : la diode conduit, V=0 V. 

dO/dt + Q/T =1I (il y a déstockage). 

1<0 et Q =0 : la diode se comporte comme une capacité (de diffusion) 
Ct. On a I=Ct.dV/dt. 


Ce modèle est justifié par l’étude expérimentale de la communication 
d'une diode. Le signal d'entrée rectangulaire est défini dans un tableau 
A(200), sur 200 points. L'intervalle de calcul DT fournit la largeur tem- 
porelle d'un point. Horizontalement, la fonction appliquée au circuit dure 
ainsi un temps 200 DT. 

l'est possible de tenir compte de la variation de l'épaisseur de la zone 
de transition en fonction de la tension inverse aux bornes de la diode 
par la relation C1 =Cm.(0,6) + V)-0:5. Ici, nous avons une capacité qui 
varie de 10 pF (pour V=0 V) à environ 2 pF (pour V=—12 V). 


LE PROGRAMME 


Ligne 70 : initialisation des calculs. 
Ligne 110 : définition du signal d'entrée. 
Lignes 140 à 230 : calculs et dessin des courbes | et V. 


Sous-programme 250 : dessin des axes. 
Sous-programme 320 : exemple de variation de la capacité de transi- 


tion. 
LES VARIABLES 
R : résistance du circuit (= 10 kQ). 
T : durée de vie des porteurs minoritaires (= 10:65). 
DT. : intervalle de calcul (10 ns). 


: capacité de transition (Ct) de la diode en inverse. 
A(200): contient le signal appliqué du circuit (en V). 
E : valeur instantanée de la tension appliquée au circuit (en V). 


V : tension aux bornes de la diode. 

| : courant dans la diode. 

(e) : charge stockée. 

DV : variation de V pendant l'intervalle de temps DT. 
DQ : variation de Q pendant l'intervalle de temps DT. 
N : position horizontale du curseur (0<N=< 200). 








——— -- me RE 
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Pour le régime rectangulaire, remplacer la ligne 110 par : 


110 : FOR N=1 TO 20:A(N) =10:NEXT 
115 : FOR N=21 TO 200:AI(N) = — 10:NEXT 


110 : FOR N=1 TO 100:A(N) = 10:NEXT 
115 : FOR N=101 TO 200:A(N) = — 10:NEXT 


La capacité dépend de la tension inverse. 


C(V) est donnée ligne 320 


RÉGIME SINUSOÏDAL 


10 REM HHHRRRE RER RER ARR EEE RER 


2C REM *# # 
JO REM # DIODE EN COMMUTATION « 
40 REM * + 


SC REM HARMAN EEE EE 
ED INK O.C:INK 1,Z6:PAPER D:PER 1 






70 R=1000C:T-0.00000! :DT=0.00090001 : 215-141 

PO REM EXEMPLE LE SIGNAL D'ENTRÉE 

9G DIM A(290) 

100 GÔSUE 320 

110 FOR N=C TO ZOG:A(N)=IO4CINESRBTEN/200)eNDT 


120 CLS 

190 GOSUR 250 

140 FOR N=0 Tü 202 

150 E=A{N) 

{0 I=(E-V)/R 

170 IF IC AND V£O THEN C=FNC(V):0G=0:DG=0:DV=IaDT/C:607T0 210 
190 IF IiO THEN V=0.6:DV=0:Dü=(I-0/7T)#DT 

190 IF IEC AND QÿC THEN V=G:0ÿ=0:DQ=(1-Q/TY#4DT 

205 IF I£0 AND Q<=0 THEN C=FNCIV) 2 G=0:DG=C:DV=I#0T/C 

216 Q=Q+D0:V=V+DV 

220 PLOT Z#N,100+4000C*]:PLOT Z4N, 300+10#V 

240 NEXT N 

240 LOCATE i,1:END 

250 REM AXES 

260 PLOT G, 105: DRAMR 450,0:DRAUR -5,5: 
370 PLOT O,J00:DRAUR 450, O:DRANE -5,5:5 
290 LOCATE 30,19:PRINT ‘I" 

J00 LOCATE JC, 7:PRIAT “4° 

315 RETURN 

J2C DEF FNC(V)=7,7E-12/{(0.8-4)49.5i 
330 RETURN 
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REDRESSEMENT 


Ce programme propose l'étude de circuits classiques destinés au redres- 
sement. Pour des alimentations conçues pour de faibles débits, on uti- 
lise un montage simple alternance (figure 1) ; sinon, on effectue un 
redressement double alternance (figure 2). On s'intéresse en général 
à la qualité de la tension continue obtenue en sortie et au courant tra- 
versant la diode. 
Le programme de base est relatif au montage de la figure 1. La tension 
d'entrée VE est définie dans un tableau A(200) ; sa fréquence est celle 
du secteur : 50 Hz. La présence de la diode D impose deux types de 
fonctionnement : 
— VE>XVS : la diode est conductrice ; ou a VE=VS. 

et1= VS + cv 

R dt 

— VE<VS : la diode est bloquée ; le condensateur C se décharge dans 

la résistance R ; ou a 1=0. 
VS __ C dVs 


s dt 





et donc 








Les différentes valeurs de VS (il y en a 200 au total) sont mises en 
mémoire dans un tableau B(200). Celui-ci permet, après l'exécution des 
courbes par le programme, de calculer, puis d'afficher, les grandeurs 
caractéristiques de VS : valeur moyenne et ondulation crête à crête. 
Le programme se modifie aisément pour permettre l'étude du courant 
let du redressement double alternance. Tous les résultats obtenus sont 
“classiques”, bien qu'il n'existe pas de solution analytique pour ce pro- 
blème. Cependant, l'intérêt d'une telle étude est qu'elle permettrait, 
par exemple, de tenir facilement compte d'une charge R variable avec 
la tension appliquée VS. On pourrait aussi étudier l'influence des imper- 
fections de la diode (résistance directe, charge stockée). 


LE PROGRAMME 

Lignes 10 à 20  :initialiation du calcul : choix deR et C (on pourrait 
mettre des INPUT). L'intervalle de calcul est pris 
égal au 1/100: de la période secteur (1/50e de 
seconde). 

Ligne 50 : définition de la tension d'entrée [enlever ABS() 


pour l'étude du redressement simple alternance]. 
Lignes 60 à 140  : calcul et affichage de VS et de | (dans la diode). 
Lignes 150 à 200 : évaluation des valeurs moyennes et de l’ondula- 
tion crête à crête. 
Lignes 210 à 240 : inscriptions sur l'écran. 


Il y a deux ‘’passages”” sur l'écran, le premier étant le régime transi- 
toire, le second le régime permanent. 


LES VARIABLES 


DT : intervalle de calcul (1/100e de 20 ms=période secteur). 

R : résistance du circuit (en ohms). 

C : capacité du condensateur (en F). 

T : RC constante de temps du circuit. 

A(200)  : tableau contenant les 201 valeurs de VE. 

B(200) : tableau contenant les 201 valeurs de VS. 

UÜ : variable de boucle locale. 

DV : variation de VS pendant le temps DT. 

MA, MI  : variables contenant, à la fin de la boucle 150-200, les 
valeurs maximum et minimum de VS. 

S : somme de 100 dernières valeurs de VS, permettant, ligne 


220, le calcul de la moyenne de VS. 





Figure 1 
Redressement simple alternance. 


1 REM HHRHHERRRRENEEEE 

2 REM # L 

3 REM # ETUDE DU + 

4 REM * REDRESSEMENT # 

5 REM # # 

b REM KEINE 

10 DT=1/(100#50) :R=100:C=0. 0005 : T=R#C 

20 DIM A(Z40):DIM B(240}:N=0 

30 CLS 

50 FOR U=0 TO 200:A(U)=ARS(GO4SIN(Z42T#)/100)):NEXT 
60 FOR U=0 T6 205 

70 D=A(U+1)-A(U) :DV=-VS#DT/T:VE=A(U) : EEU)=VS 

80 IF ARS(D)>=ARS(DV) AND 130 THEN VS=A(U: ELSE VO=VS+IN 
100 


90 IF VE:VS THEN VO=A(U) 

O0 I={VS/RI+(D4C/ET) 

IF I£C GR VSSVE THEN I=0 

PLOT 204244, 200+VE:PLOT 20H20, ZODHVE 
PLOT 204240, 304104] 

NEXT UsNEN+i:IF NE THEN CLS:G0T9 &û 
REM CALCULS 

-MA=9G:NI=90 : 520 

FÜR U=0 TÜ 206 

180 IF B(UÏ£MI THEN MI=E(U) 

190 5=5+B{i) 

206 NEXT U 

220 LOCATE 30, SPRINT VOS": PRINT USING'A#4 HS: CINT(S/200)3: 

235 LOCATE 30, 2:PRINTDV=" 3: FRIAT USING AH. #5: (INTUMA-MI) je: PRINTS UE 


CO = OO Un + Ci PO pa 
5 1 À © © © © Q 








240 LOCATE 40,3:END 





Figure 2 
Redressement double alternance 
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REDRESSEMENT SIMPLE ALTERNANCE 
R=100 

C=500 uF 

(enlever la fonction ABS, ligne 50). 











REDRESSEMENT DOUBLE ALTERNANCE 
De haut en bas : 

vs 

VE 


EE 





DETECTION SYNCHRONE 


Ce programme est à comparer au précédent. Il permet de démoduler 
une onde multipliée en amplitude. Le circuit de principe est le suivant : 


MULTIPLIEUR 
VE = ———- x 
Vu ti + sin wgt) sin wyt. 
HF Locale 


SIN (ot + w) 





On démontre que : 


Vs = Vu cos p+ — VM cos @ sin wBt — 
+ -VM(sin wBt +1) cos (2uHt +4) 


Le détecteur se compose donc simplement d'un multiplieur (que l'on 
appelle mélangeur équilibré en HF) ; sa sortie fournit une composante 
continue, l'onde BF démodulée et une nouvelle onde modulé en ampli- 
tude, de pulsation 2wH que l'on élimine à l’aide du filtre RC (passe bas). 
‘On comprendra aisément que ce type de démodulation est beaucoup 
plus linéaire que le précédent. Rappelons que le rapport signal/bruit est 
aussi amélioré de 6 dB. L'onde sin (wHt+) doit être absolument 
synchrone de la porteuse. En effet, leur déphasage + intervient direc- 
tement dans l'expression de l'onde démodulée, et il est hors de ques- 
tion qu'un tel déphasage varie en fonction du temps. A titre d’exem- 
ple, on peut voir sur les photos l'influence de la valeur de # sur l'aspect 
de l’onde démodulée. 

Le programme est simple. L'expression de l'onde modulée VE est éta- 
blie, ainsi que celle de VS. Le circuit RC est mis en équation de façon 
classique. 
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VS=RI+V 


SALUE 
dt 
d'où av= VS Vidt 
RC 


Il Y a donc deux ‘‘passages”” sur l'écran : seul le second a été photo- 
graphié, le premier, en effet, est affecté du régime transitoire du 
système. 


LE PROGRAMME 
Lignes 60 à 120 : définition des paramètres. 
Lignes 160 à 240 : calcul des différentes tensions. 


Sous-programme 270 : tracé des axes et inscriptions sur l'écran. 


LES VARIABLES 


BF _: fréquence de l'onde basse fréquence (en Hz). 

WB : sa pulsation. 

HF _: fréquence de l’onde porteuse (en Hz). 

WH : sa pulsation. 

PH __: déphasage (en radian) entre la porteuse et la HF locale. 
M: taux de modulation (0<kM<1). 

DT : intervalle de calcul (en s). 

VM _: amplitude de l'onde HF. 

FC __: fréquence de coupure adoptée pour le filtre passe bas. 
RC __ : valeur correspondante de la constante de temps. 

T : variable temps. 

L : compteur du nombre de ‘passages’ sur l'écran. 

DV  : variation de V pendant DT. 

VB  : expression de l'enveloppe de l'onde modulée. 

VE _: expression de l'onde modulée. 

VS  : tension de sortie du multiplieur. 

V __: tension de sortie du filtre passe bas. 


Les courbes présentées montrent l’aspect de la tension BF (en haut) 
et de la tension démodulée (en bas) en fonction du déphasage entre 
porteuse et oscillateur local. On constate que la tension de sortie est 
la somme : 

— d'une tension continue, 

— de la BF modulée, 

— d’une onde modulée en amplitude, de pulsation 2 wHF. 

Il n'y a pas détection pour PH=902. 
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Il 
fi 
wi 
Lu 


fu 
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————. FÉÉACE-AC 


_— —. FE rs ini. ES 2% 









































— LIST — 
mi re = BSITT ASS s— LE 





10 REM HAS REH RENÉE REREERERSE 


20 REM +# * 
JD REM * DETECTION SYNCHRONE « 
40 REM * # 


50 REN FF#EIARRHIIAREEFERESRERTE 

ëÿ LS 

70 EF=4000:KP=S#PI#Er 

EC HF={O0OCC:HH=SAF TIRE 

90 M=t:D1=0. 000001 :VH=i 

100 T=G:FC=0COC: INPUT PRSESCZ(EN DEGRES)" 8: PH=049:/100 
110 RC=j/(2#PIxFC; 

120 CLS 

149 GOSUR 270 


Q LaL41 


319 
239 








G STDT 

ü VB= VE € LES IN CNBAT)) 

C VE=UMRSINEMHET IEC LEMESININSET)) 
200 VE=VERSINCNHET+FH) 

EC DR Lier 

220 PLOT ZaÙ, 195+00eVr 

ED PLOT za, IGH2OCGRV 

240 NEÏT ü 
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Z THEN LGCATE 1,13:END 

260 6070 {20 

270 REM AŸES 

280 FOR N=0 TO 240 STE0 35 

E9C PLOT Z4N, {ic 

300 NEXT 

330 LOCATE 20,1:PRINT "PHASEZ* SUGING SAYS {Bies PRINT * DE 
350 RETURN 





DETECTION 


Ce programme, presque identique au précédent, s'attaque à un pro- 
blème pourtant différent : la détection d'une onde modulée en ampli- 
tude, c'est-à-dire la démodulation. Il permet d'obtenir et d'étudier la 
forme du signal recueilli en sortie du détecteur. Il met en évidence le 
fait que la démodulation ne peut pas être parfaite quand le taux de 
modulation est de 1. 

Le circuit détecteur étudié est identique à celui de la figure 1 du para- 
graphe ‘‘’Redressement'”. L'expression d’une onde modulée en ampli- 
tude est donnée par : 

VE : VM Sin wHt x (1+m Sin wBt), dans laquelle : 

— VE est la valeur instantanée de l'onde modulée. 

— VM l'amplitude de l'onde (sinusoïdale) non modulée, c'est-à-dire 
quand m=0. 

— «wH est la pulsation de la porteuse (HF). 

— wB est la pulsation de la tension de modulation (BF). 

— m est le taux de modulation (o<m=<1). 

Le condensateur C se charge à la valeur crête de VE pendant ses maxi- 
mums, puis se décharge dans la résistance R pendant les ‘’Creux'’ de 
la modulation. La tension VS est donc une image de l'enveloppe de 
l'onde modulée ; on recueille ainsi en sortie du détecteur une tension 
basse fréquence (BF). Celle-ci n'est pas parfaite pour m=1 (cas du pro- 
gramme). (On démontre en effet qu'il faudrait avoir : 


1 m : 
> ce qui est impossible pour m=1). 
TOR EE à rer 2 Là LS 


Le programme est difficilement exploitable sur un moniteur vidéo, quelle 
que soit la qualité de la haute résolution du micro-ordinateur employé. 
En effet, pour que le phénomène étudié soit clairement visible, il fau- 
drait que l'axe OX ait une définition supérieure à 1000 points. L'’impri- 
mante convient ici particulièrement, puisque son axe OX est infini. 
Ce programme est très court. 
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Lignes 60 à 90  : définition des constantes. 
Lignes 100 à 180 : boucle de caicui. Le système est en équation lignes 
130, 140, 150. 


LES VARIABLES 


BF _: fréquence de l'onde basse fréquence (en Hz) ; ici BF =2 kHz. 
WB_ : pulsation de l'onde basse fréquence (en rad/s). 
HF __: fréquence de l'onde haute fréquence (en Hz) ; ici HF = 50 kHz. 
WH_ : pulsation de l'onde haute fréquence (en rad/s). 


M : taux de modulation ; ici M=1. 

RC: constante de temps du circuit détecteur (en s) ; ici RS = 100 us. 
DT  : intervalle de calcul en s ; ici DT =0,5 us. 

VM__: amplitude de l'onde HF non modulée. 

T : variable temps. 

VE : valeur instantanée de l'onde modulée. 


VS  : tension de sortie du circuit détecteur. 


40 R M Li # 

HRFRRREPHEERE 

F= 2000: HE=CEPIARE 

BC HF=SO0CO:HH=Z4PI#HF 

96 Mi :RC=S, 0001 : 0720. COCOCOS : VM=1 

106 T=Ù 

110 CLS:INK OO:INK L,ZH: PAPER CiFEN 
115 PLOT GC,ZCC:DRANS 64C,C 

120 FGR N=0 Tü 640 STEP 0.5 

130 VE=VMESINCMHET IEC LEMESINUMBXT) ; 
140 IF VES=VS THEN DUSVE-VG ELGE Tiÿ=-VS4nT/ 20 
150 VS=VS +0 

160 T=T#DT 

176 FLOT N,200+90#VS 

180 PIGT N,200+90HVE 

196 NEXT K 

200 FOR 2=1 TO SOOCG:NEXT 1 

216 CLS:6670 129 


ss tn 
Co JE mue ] 
td 2 


110 


“ 


LORS 








COR 





Chapitre VI 


ABAQUE DE SMITH 


111 


Cet abaque, imaginé par l'ingénieur Smith, permet de résoudre un grand 
nombre de problèmes relatifs à la transmission de signaux HF. 

Il s’agit d’un cercle, à l’intérieur duquel on peut représenter n'importe 
quelle impédance, un déplacement le long d’un câble de transmission 
se ramenant à une rotation dans ce plan. Le programme que nous pro- 
posons utilise la méthode classsique d'exploitation de l’abaque, mais 
par l'emploi de ses propriétés analytiques et non graphiquement. 
Cette abaque est un cercle, dont tout point intérieur permet de repré- 
senter une impédance normalisée. Soit Z l’impédance étudiée et Zc l'im- 
pédance caractéristique de la ligne employée ; 


Z : : , ; 
alors Fe un nombre sans dimensions, s'exprimant par : 
C 


+ =R+jX, R représentant la partie résistive, X(> <O) la partie 


réactive. 


— un R donné correspond à un cercle dont le centre est sur l'axe des 
x d'abscisse. 


et d’ordonnée O. 





R 
R+1 





. 
Le rayon de ce cercle est 

R+1 

— un X donné correspond à un cercle dont le centre a pour abscisse 
1, pour ordonnée 


_. ; le rayon de ce cercle est 1/X. 


Ceci permet de tracer la figure 1. 

D'autre part, le grand cercle est gradué en longueur d'onde, un tour 
complet correspond à )/2. 

Rappelons la manière dont on utilise une telle abaque. 

Soit un générateur attaquant une ligne d'impédance caractéristique ZC, 
fermée au bout d’une longueur | sur une impédance ZR. Quelle est l'im- 
pédance Zg que ‘‘voit’’ le générateur ? 
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La méthode consiste à : 
— placer sur l’abaque le point représentatif de : 


——=R+jX (point A). 


(Sur le dessin, nous avons ici R=0,33 ; X= +0,5). 


— Faire tourner ce point autour du centre O de la figure, en lui faisant 
décrire un arc « tel que : 


4rl 
Q= 


À 


— Placer le point B ainsi obtenu et “‘lire”’ ses caractéristiques (ici, on 
a environ : R'=2, X'= +1,5). 


— En déduire ZG=(R"+jX")Zc. 


Cette méthode graphique est extrêmement pratique et permet d’ailleurs 
de résoudre une infinité d'autres problèmes. En fait, mis au point à une 
époque où les calculatrices n’existaient pas, c'était le seul moyen rapide 
d'obtenir la solution. 

Maintenant, nous pouvons reprendre les équations: des cercles et trai- 
ter le problème de façon purement numérique. Ncus en faisons grâce 
au lecteur et lui épargnons la démonstration. Les équations que nous 
fournissons apparaissent dans tous les traités un peu âgés sur l'aba- 
que de Smith. 


— Connaissant R et X, on peut en déduire la position du point corres- 
pondant dans le plan. On a : 


FE 2X 
VO RE 
Position du point A : 
Yo(R+1) 
X0O= 1 - ——— 
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— Une rotation d'angle @ dans ce plan amène aux nouvelles coordon- 
nées : 


x1=x0 Cos 8—yo Sin 8 
y1=x0 Sin 8 +yo Cos 8 


— Une transformation inverse conduit aux résultats : 


x' 2y1 
x1?+y12—-2x1+1 


x'(1—x1) 
y1 


R'= 1 





Légende de la photo 
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ABAQUE DE SMITH 
IMPEDANCE DU CABLE:7 75 
RESISTANCE DE LA CHARGE:?7 
REACTANCE DE Lh CHARGE : ? 166 
LONGUEUR DE LA LIGNE CEN LONGUEUR 
D'ONDE):7 .45 


1G0 





RG=-49.3 OHMS 
XG=67.7 CHMS 


ART ne 


or : 
DES 5 se 
_ : : 
2* à LES 
Par È A 
. es PRES 





‘ : ; 
s f ra, Se 
F ‘ - “. — 
J j Fe ES #-- 3 
1 à 
: ‘ | - Rue 
H Î cb 
b à 1, . : ee ‘ 
a À ; À, Ne 
\ ; Se LS ee 7 
NN ai” Mens + 
ne É— 
DER est l 





Exemple d'exécution du programme 1. 
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PROGRAMME 1 


10 REM #RHHRHRRERREREREE 


20 REM # 4 
30 REM # ABAQUE DE SMITH # 
40 REM * Le 


50 REM HEAR 
70 GOSUR SOC 

80 RL=R:XL=X 

95 CLS 

100 FOR N=1 TO 6 

110 READR 

120 FOR J=0 TO Sé0:1=J#PI/180 


130 PLOT 200+(1804R/(R+1))+(1804C0S(1)/{R+1)),1GC+(1BCKSIN(I)/{R41)) 


140 NEXT J:NEXT N 
150 DATA 0,0.33,0.66,1,2,4 
160 FOR N=i TO 3 
170 READ X, 70 
180 T=70 
190 PLOT 380-180C0S(T)/X, 180-(180/X)+1B04SIN(T)/X 
200 T=T4P1/128 
210 IF T>=PI/2 THEN 220 ELSE 190 
220 DATA 2,-.638136, 1,0.024543,0,5,0. 643592 
230 NEXT N 
240 FOR N=1 TO 3 
250 READ X, T0 
260 T=-PI/2 
270 PLÔT 380-180#C0S1T)/X, 180+(180/X)+1B0RSINET)/X 
200 T=T4PI/128 
290 IF T>=T0 THEN GOTO 300 ELSE 270 
300 NEXT N 
310 DATA 2,0.662679,1,0,0.5, -0. 443592 
320 REM CALCULS 
330 R=RL:X=XL 
340 R=R/2C:X=X/2C 
350 Y=2#X/((XA234(R41)A2) 
360 X=1-Vx{R#1)/X 
380 LOCATE (196+180#X)/16,25-((172+180#V)/16): PRINT" +" 
400 XN=X#COS(L)-Y#SIN(L) 
410 YN=XASINIL)+Y#COSIL) 
20 NX=ZEYN/ CXNA2EYNAD4 1 (ZHYN) ) 
430 NR=(NXKE1-AN)/YN)-1 
440 LOCATE (196+180#XN)/16,25-((172+1RC#VN)/16) :PRINT"a" 
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62 USING He NREZCS PRINT 


460 LOÜCATE 35, 
5 





470 LGOCATE 45, NT * 

490 GÜSUB 9100 

490 GOTQ 490 

500 REM INTRODUCTION DES DONNEES 

5i0 CLS 

320 PRINT: PRINT:FRINT * AEAQUE DE SMITH 


3930 PRINT: INPUT"IMPEDANCE DU CAFLE:":2C 

540 PRINT: INPUT'RESISTANCE DE LA CHARGE: IR 

330 PRINT: INPUT'REACTANCE LE LA CHARGE: "IX 

560 PRINT: INPUT'LONGUEUR DE LA LIGNE (EN LONGUEUR 
370 L=-4*PI#L 

580 RETURN 


Exemple d'exécution du programme 2. 


ABARUE DE ZMITH 


D'ONDE sel 


LONGUEUR DE LA LIGNE: L=9.H9S LONGUEUR D’ 
ONDOE 

Ze REDUITE Ze KEDUITE 
5.5 #+ 3 1.14 + 1.S2 
G.6 + 1.34 #+  J L:ST 
6.7 + J 1.S2 + y 1.44 
6.3 + J 32:59 DE | L. SX 
GG. 3 Fe -ÿ 1.35 + 1.16 
L + J 1.936 + J 1.51 
Lu + à 2.22 + J 8.285 
L:Ee UE 2.18 Ÿ 0.7 
1.5 + 2.26% + .j 6.352 
1.4 + J E -«33 + J 6.35% 
1.5 + J 2.39 + J 9.23 
1.6 LAURE | 2.423 + J 5.08 
1:27 + 2.47 — J 9.25 
L:6 LE 2.493 =) 0 2.19 
L:9 + .j BR: Si = 3 6.233 
2 + J 2.S1 — J 6.46 


PROGRAMME 2 


10 REM HAHHERRHIERENS ER 


20 REM * * 
JO REM # ABAQUE DE SMITH # 
40 REM *# + 


50 REM ###H HN RNIERERIREES 
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Chapitre VII 


PLL : PRINCIPE 
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CONSTITUTION D'UNE BOUCLE 
A VERROUILLAGE DE PHASE 


Une boucle se compose de quatre parties : 

— un comparateur de phase (CP), chargé d'élaborer la différence entre 
les phases de deux tensions (figure 1). 
Suivant les cas, la caractéristique de transfert peut être triangulaire 
(figure 2) ou sinusoïdale (figure 3). 
En régime linéaire, nous écrirons dans les deux cas : 


Vdd/2 
x/2 
En régime non linéaire, la caractéristique figure 2, discontinue, ne se 


prête que difficilement à une exploitation mathématique ; nous adop- 
terons, dans ce cas, l'équation de la caractéristique figure 3, soit : 


KPp=V1/p= = Vdd/7x 


V1 =Vad/2 sin (@i— pr) = Vadd/2 sin +. 


— Un filtre passe-bas chargé de : 
— définir la zone de capture Fc (voir plus loin), 
— assurer la stabilité du système, 
— filtrer la composante à wi+wr, générée en même temps que le 
signal utile dans le CP. 


— Un VCO, oscillateur commandé en tension, délivrant une tension de 
sortie, de fréquence proportionnelle à sa tension d'entrée. Sa trans- 
mittance est appelée Ko. 

La définition des différents termes est donnée figure 4. 


— Un intégrateur, ‘’immatériel’”, dans la mesure où il n’a aucune pré- 
sence physique, signalé dans le schéma figure 4. En effet, la gran- 
deur de sortie du VCO est la pulsation d'un signal dont on utilise 
à l'entrée de la boucle la phase, c’est-à-dire l'intégrale par rapport 
au temps de la pulsation : pr= wdt. 


FONCTIONNEMENT DE LA BOUCLE 


1°) Régime transitoire : cas d’un échelon de fréquence. 

En l'absence de signal d'entrée, la boucle délivre sa pulsation centrale, 
soit wO. À l'apparition du signal, tel que sa pulsation soit comprise entre 
&0 — wL et w0 +w, il naît, à la sortie du comparateur une tension, varia- 
ble, fonction de la différence des phases des signaux d'entrée et de 
réaction, de la forme : v1=KP (@i— or). 

Cette fonction croît dans le temps si wi>wO (en négligeant bien sûr 
toute non-linéarité pour le moment). En l'absence de filtre passe-bas, 
ona:v2=vi(>0)et donc le VCO voit sa pulsation de sortie augmen- 
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ter et passer à wNO. Ceci se traduit par une moindre augmentation de 
v1(=KP.[wit - wNOt]}), mais une nouvelle croissance de wNo. Ce phéno- 
mène se poursuit jusqu'à la capture, moment précis où wNO =wi ; à cet 
instant là, v1 est devenu constant, et correspond à un certain dépha- 
sage (constant) entre vi et vr. La boucle est verrouillée. 

En pratique, ce régime transitoire n’est pas aussi simple que cela, pour 
plusieurs raisons. 


2°) Les imperfections d’une boucle réelle. 


— |l y a d’abord un terme parasite à la sortie du comparateur de phase 
dû à son principe de fonctionnement même. En effet, que celui-ci 
soit constitué d’un OÙ exclusif ou bien d’un multiplieur différentiel, 
le signal d'erreur est constitué du produit des deux tensions vi et vr: 


vi= Vi sin wit 
vr=Vr sin wrt 


vives SE [Sin (wit — wort + … — COS (wit +wrt)] 


Nous voyons donc apparaître un terme en wr+wi, en plus du terme 
utile en wi—wr. || se superpose donc au signal d'erreur v1. Cepen- 
dant, en régime permanent, il ne modifie pas la pureté spectrale de 
v1, dans la mesure où sa valeur moyenne est nulle sur une demi- 
période de vr. 

Ensuite, le comparateur de phase n’a pas la caractéristique rectili- 
gne indiquée figure 4. Au-delà de x/2, il y a en effet ‘‘retournement"" 
de la caractéristique v1(9), que celle-ci soit linéaire (OU exclusif) ou 
sinusoïdale : la réaction devient alors positive jusqu’à 3&/2, néga- 
tive jusqu'à 57/2, etc. Le système n’a cependant pas de point de 
fonctionnement stable dans les zones à réaction positive : il s’en 
éjecte de lui-même, pour se retrouver dans une zone à réaction néga- 
tive contiguë. 


— Enfin, la dernière difficulté provient du fait que, lors d’un échelon 
de fréquence, par exemple, on ne maîtrise que rarement la phase 
initiale avec laquelle se présente le signal à t =0o. Dans ce cas, seule 
une étude statistique permet de prévoir les limites des lieux de vi(t) 
et de connaître ainsi les valeurs extrêmes que peuvent prendre les 
paramètres de fonctionnement, tels que temps de réponse, etc. 


3°) Le filtre passe-bas. 

Bien que la réponse de la boucle soit satisfaisante en l'absence de fil- 
tre, il est presque nécessaire d’en inclure un dans la boucle et ceci pour 
plusieurs raisons. 

Remarquons d’abord que, contrairement à beaucoup de systèmes asser- 
vis, on ne s'intéresse ici que peu au signal de sortie vr. On utilise en 
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fait les tensions V1 ou V2 qui sont les images, en régime établi, de la 
fréquence du signal d'entrée Vi. Or, sans filtre, on a pu constater pécé- 
demment que, au signal utile, se superpose une ondulation de pulsa- 
tion wi+wr. Celle-ci perturbe le régime transitoire et doit de toutes 
façons être éliminée de V1. Elle constituerait un ‘‘bruit”” inadmissible 
en démodulation de fréquence. L'utilisation d’un filtre a donc pour but 
de faire disparâitre ce signal gênant et, par là même, tout signal inter- 
férant parasite. Cependant, ce filtre ne doit pas pour autant faire dis- 
paraître l'erreur ‘‘utile’” en wi —ur, responsable du verrouillage de la bou- 
cle : d’ailleurs, son atténuation est responsable de l'existence d’un nou- 
veau paramètre, appelé fréquence de capture (fc). Pour qu'il y ait ver- 
rouillage, il faut que la fréquence incidente (fi) soit comprise dans un 
domaine de largeur 2fc entourant fo. Cette contrainte disparaît dès qu'il 
y a eu verrouillage. Le choix de la constante du filtre agit donc directe- 
ment sur fc. 

Enfin, le fait d'introduire un tel circuit dans la boucle la transforme en 
un système du second ordre : il faudra alors prendre garde à la stabi- 
lité. 


BOUCLE SANS FILTRE PASSE-BAS 


Le schéma est alors très simple (voir figure 1). En boucle fermée, le 
système est du premier ordre et donc toujours stable. 


1°) Régime linéaire. 

Soit Vad/r la transmittance du comparateur de phase. 
Soit 2wL/Vad la transmittance du VCO. 

Le gain de boucle s'écrit : 











did Vad  2ét else 20 à 2h 
7 Vdd p T p 
un 2owL/T 
AIRE 1+pr/2ul 


La constante de temps de ce système est 7=7/2wL et le temps de 
réponse à 5 % vaut : 


3T 
2wL 





37 = 
On peut éliminer la fonction SIN de la ligne 220 pour étudier la boucle 
dans cette hypothèse. 


2°) Régime non linéaire. 
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C'est naturellement le fonctionnement normal de la boucle : la résolu- 
tion numérique présente tout son intérêt. Le programme permet en effet 


de 


tracer la réponse à un échelon de fréquence, en tenant compte de 


la phase initiale aléatoire avec laquelle se présente l'entrée vi à t=0. 
Il y a plusieurs possibilités d’'détude : 


LE 


Lig 
Lig 
Lig 
Lig 
Lig 


Tracer la réponse pour plusieurs valeurs fixées de cette phase (PHI 
en degrés) : 0°, 30°... Comme prévu, tant que PH est compris entre 
— 90° et +902, le système démarre d'une zone à réaction néga- 
tive et à un comportement linéaire ; en dehors de cet intervalle, la 
réaction est d’abord positive ; cela se traduit, au départ, par un écart 
à la valeur finale grandissant dans le temps, puis passage de V1 à 
la valeur — 1, enfin retour dans une zone à réaction négative. Le phé- 
nomène est spectaculaire pour PHI=100° ! 


Tracer la réponse pour des valeurs aléatoires de PHI et superposer 
sur l'écran toutes les trajectoires correspondantes. Cela permet ainsi 
de mesurer, dans les conditions les plus favorables, le temps de 
réponse de la boucle. Il faut alors remplacer la ligne 100 par : 


LET PHI= — 360 *RND + 180 


Ilest possible, enfin de calculer par programme le temps de réponse 
de la boucle pour de valeurs aléatoires de PHI, puis de tracer un dia- 
gramme à barres permettant de connaître le temps de réponse le 
plus probable. 


PROGRAMME 


nes 60 à 120 :initialisation' des paramètres. 

nes 140 à 180 : tracé des axes. 

ne 190 : introduction de la valeur initiale de la phase. 

nes 200 à 250 : calculs. 

nes 260 à 330 : introduction de la valeur désirée pour la phase ini- 
tiale. La nouvelle courbe se superpose alors à celle 
obtenue pour une phase initiale nulle. 


LES VARIABLES 


: fréquence centrale du VCO. 

: fréquence de verrouillage. 

: fréquence d'entrée. 

: tension initiale à l'entrée du VCO. 
: phase initiale à la sortie du CP. 

: fréquence instantanée du VCO. 

: intervalle de calcul. 
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S : valeur de PH désirée à t=0 (introduire 4 chiffres ; l’exprimer en 
degrés). 


Comportement de la boucle du premier ordre, suivant la valeur initiale 
de la phase. 


REMARQUES SUR LE PROGRAMME 


1°) Comme on a pu le constater, nous avons éliminé, par souci de sim- 

plification, l’état de repos de la boucle. En pratique, un comparateur 

de phase est constitué par un système qui délivre au repos (fréquence 

d'entrée Fi=fréquence centrale du VCO=F0, ou bien pas de signal 

d'entrée) une tension de valeur moyenne égale à la moitié de la tension 

d'alimentation Vdd. Ceci correspond à un déphasage de +x/2 entre les 

deux entrées du CP. Nous avons donc posé, à l'équilibre, V1=V2=0 

(au lieu de Vad/2). 

#=0 (au lieu de x/2). Les notations utilisées sont les mêmes que dans 

le texte sauf : 

— PHI= phase initiale en degrés. 

— FNO=fréquence instantanée du VCO. 

— DT=intervalle de calcul, pris égal à 1/8 de la période 1/Fj du signal 
d'entrée. 


2°) L'exemple proposé est le suivant : 

FO = 50 kHz 

FL= 20 kHz 

F;i=69 kHz (valeur quelconque entre 50 — 20 = 30 kHz et 50 +20= 
70 kHz pour qu'il y ait verrouillage). 


3°) On constate que les résultats ‘‘linéaires’’ sont à peu près applica- 
bles lorsque la phase initiale est comprise entre — 90° et + 90°. Dans 
ce cas, en effet, le système part d'une zone où il y a effectivement 
contre-réaction (AVI/Ag > 0) et donc n’a aucune raison d’en sortir. Par 
contre, si —270<p< —-90 et +90<gp< +270°, le système voit son 
gain dynamique changer de signe (AV1/A6 <0) : il y a réaction posi- 
tive. Sur l'écran, on constate que le point de fonctionnement s'éloigne 
d'abord de l'équilibre final (c'est-à-dire Vi=0,95 dans l'exemple numé- 
rique choisi), passe par Vi= —1 (soit : y =270°), puis retourne au régime 
linéaire, avec une phase initiale à ce moment-là de +270°. Dans de 
tels cas, le temps de verrouillage peut atteindre 8 cycles de Fi. 


4°) Nous pouvons également remarquer que, même en régime de réac- 
tion négative, le temps de verrouillage de la boucle est plus important 
que ne le prévoit la théorie de son fonctionnement linéaire. En effet, 
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dans ce cas, et à cause de la caractéristique sinusoïdale du CP, le gain 
dynamique de la boucle diminue et tend vers 0 lorsque Fi se rappro- 
che de Fo=FL (ou Fo —FL) car, alors la limite finale de Vi est +1 (ou-—l), 
valeur pour laquelle la pente de la fonction sinus tend vers 0. 


10 REM RAREREIREREREE 


20 REM + € 
JG REM # PPL D'ORDRE 1 + 
40 REM + * 


50 REM HÉROS 

60 FO-SCOOG:REM FREQUENCE D'ENTRÉE 

70 FL=200C0:REM FREQUENCE DE VERROUILLASE 

80 FI=69000:REM FREQUENCE D'ENTREE 

90 VI=O:REM TENSION INITIALE 

100 PH=C:REM PHASE INITIALE EN DEGRES) 

110 FX=FD 

120 DT=1/(424F 1) 

130 CLS 

140 PLOT 19,190:DRAWR 400,0 RAENE ee 1DRAKR 5,-S:CRAKR - 
150 PLOT 10,0:DRAWR O,J80:DRANR -5,-S:DRAUR 5 | SsDRAUE : 
160 FOR N=10 TO Z10 STEP 3 

170 PLOT ZaN,10:PLOT 24N,270 

180 NEXT 

190 PH=PH*PI/180 

200 FGR I=0 T0 200 

210 PH= PRRSPETAIEl -FX) 

220 VI=SINIPH 

230 Sr 

240 PLOT 10+2+1,100+180#V1 

250 NEXT 

ZbO LOCATE 29,1:INPUT *FHAGEZ?:P 

270 LOCATE 29,1:PRINT" . 

290 LOCATE 25,20:PRINT *PHASE=":P:°Deg.* 

300 T=INT{200+07T#{ 000000) 

310 LOCATE 25,18:PRINT "ECRAN=":T:"E-9 
30 PH=P:FX=FO:GO0TO 150 


LA 
4 


Un Lun 
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Figure 1 : En posant : 
vi= Vi sin (wit + çi) 

vi= Vi sin @i 

et Vr= Vr sin ort 
Vr=Vr sin @r 

on a V1=Kp (@i— or). 





Figure 3 : Comparateur de phase à caractéristique sinusoïdale. 
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Figure 4 : Schéma synoptique 
d'une boucle à verrouillage 


de phase. 





Figure 5 : Temps de réponse de 


la boucle en fonction de la 


phase initiale. 
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BOUCLE A VERROUILLAGE DE PHASE 
DU SECOND ORDRE 


Le programme précédent rappelle le principe de fonctionnement d'une 
boucle à verrouillage de phase. Ce composant joue un rôle important 
dans les systèmes de communication moderne. En particulier, il est 
aujourd'hui universellement employé pour la démodulation de fréquence 
dans les récepteurs radio. Les notations employées dans le texte sont 
les mêmes que précédemment. 

Le programme sert principalement à étudier le régime transitoire pour 
un signal d'entrée FSK. Ce type de modulation est utilisé pour trans- 
mettre des informations binaires. Un ‘‘O'" est représenté par une cer- 
taine fréquence Fo (qui sera justement la fréquence centrale — ou de 
repos — de la boucle étudiée) et un ‘’1’” pour une autre fréquence Fi. 
Nous proposons un exemple numérique, en liaison avec une manipula- 
tion qui a été effectuée au laboratoire. Cela nous permettra de con- 
fronter résultats numériques et expérimentaux. 

Il faut donc fournir successivement au programme : 


— la fréquence centrale Fo de la PLL, 
— la fréquence de verrouillage FL. 


Ces deux quantités sont définies par construction, en fonction des 


caractéristiques du VCO. 
Le choix est ensuite donné entre deux types de filtres passe-bas. 


R 
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Vs 


GC 1 
V FC. nn 
I - : 2T(R4 + Roc 


Le filtre de type À donne en général une stabilité plus faible à la boucle 
qu'un filtre de type B. 

Il faut enfin fournir la valeur de la fréquence d'entrée FI. 

Le programme dessine alors l'aspect de la tension V2 (cas A) ou V2 
et Vs (cas B). Le signal en VS est, dans ce cas, mieux filtré qu'en V2. 


LE PROGRAMME 


Il permet donc d'étudier la réponse d’une boucle, dans le cas d’un éche- 
lon de fréquence. 


Lignes 80 à 120 et 190 : entrée des données. 
130 à 180 et 200 à 210 : initialisation des variables. 
Lignes 220 à 330 : calcul et affichage des valeurs. 


Sous-programme 350 à 420 : tracé des axes ; inscriptions sur l'écran. 


LES VARIABLES 


FO : fréquence centrale. 

FL : fréquence de verrouillage. 

H : détermine le type de filtre. 

FC : fréquence de coupure du filtre passe-bas. 
VDD  : tension d'alimentation. 

RC : constante de temps du filtre. 

FXO  : valeur instantanée de la fréquence du VCO. 
v1 : tension d'entrée du filtre. 

V2 : tension de sortie du filtre (aux bornes de C). 
DT : intervalle de calcul. 

PH : valeur de la phase à la sortie du CP. 


DV2  : variation de la tension aux bornes de C pendant DT. 
Vs : tension d'entrée du VCO. 
KV : transmittance du VCO. 
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DE <CUPURE 


D'ENTREE : ? 








Données numériques valables — sauf contre-indication — pour toutes 
les photos. 





DE L’'ECRAN:G.4E-TERSEC. 





—————— ve 











us 























REGIMES TRANSITOIRES DANS UNE P.L.ZL. 
FREGUENCE CENTRALE:T7 SEA 

FRERUVENCE DE VERROUILLAGE:7 e2E4 

FILTRE RC «(12 OU FILTRE R1,R28,C €<2):7 2 


FREGQUENCE DE COUPURE C1/2#PI#RC)?7 1E3 
Re/(R1+R2)2=-7?7 6.1 
FREGUENCE D’ENTREE:?7 6HE3 
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LARGEUR DE L'’'ECRAN:9.6E-G35SEC. 


anse 




















Fi=57,8 kHz filtre (1) 





REGIMES TRANSITOIRES CANS UNE P.L.L. 
FREGUENCE CEMNTRALE:7 SE4 

FRERUENCE DE UVUERROUILLAGE:7 eE4 

FILTRE RC <C1) OU FILTRE R1,R2,C €C2>:% 1 


FREGUENCE DE COUPURE C1/2#PIKHRCD7 1EZ 


FREGUENCE D’'ENTREE:?7 S7.SE2Z 





+ LARGEUR DE L'’'ECRAN:G.SE-OD3SEC. 
{ 








F1= 60 kHz filtre (1) 
absence de verrouillage 
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FILTRE TYPE (2) 





F1= 55 kHz 

















REGIMES TRÉMSITOIREZ DANS UNE ?.L.L. 
FRERUENCE CENIPALE:?7 SE4 

FRERMENCE DE VERERQUILLAGE:?T LE 
FILTRE RC <«<Zz)> OU FILFPE FL,R2Z,C <£r:T 
FRÉQUENCE DE COUPURE (C1/2MPI#RC)T 152 
FRERUENCE D'’ENTREE:7 SSE2 

… LARGEUR DE L’'ECRAN:G.1E-GDESEC. 

Le " UE 
os : 

US 
Fi=57,8 kHz 

REGIMES TRANSITOIRES DANS UNE P.L.L. 
FREGUENCE CENTRALE:? SEA 

FREGUENCE DE VERROUILLAGE:7 SEA 
FILTRE RC (12 OU FILTRE R1,R2,C (2):9 
FREGUENCE DE COUPURE (1/2PI#RCD7 1E2 


+ 


* 


s 





—— 


Re/(R1+RE2)2=7 O.1 
FREGUENCE D'’'ENMTREE:7 


[= 
red 


7.SE 


2 
e 





LARGEUR DE L'ECRARN:S.SE-G2SEC. 


LS 
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QI 


ET 
REGIMES TRANSITCOIRES ANS UN P.L.L. 


FRERUENCE CENTRALE:7 SE4 
FREGUENCE DE UVUERROUILLAGE:7 ex 
FILTRE RC <«12 OU FILTRE R1,R2, €<2):7 2 


IÉRC:? 1E5 


% 


FRÉQUENCE DE CCUPURE (1/2 
Re/(R1+RE)=7 O.1 
FRÉQUENCE D'ENTREE:T SV.SE 


(è 


4 LARGEUR DE L'ECRAN:G.SE-O023SEC. 


U2 




















TRANSITOIRES EANT YJNE P.L.L. 
FREQUENCE CENTRALE:?7 SE4 

FREQUENCE DE VERROUILLAGE:?7 2E4 

FILTRE RC «(12 OÙU FILTRE R1,R2,C <22:7 2 


FREQUENCE DE COUPURE (C1/2%PI#RCD?9 1E5 
R2/CR1+R22>-?7 9.1 
FREQUENCE D'’ENMTREE:?7 S1E3 


+ LARGEUR DE L'’'ECRAN:OD.4E-Ge2SEC. 


= ————_—_—_—_—_—_—_—_—_—_—_——— #2 
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SEC, 


RINT 
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Oscillogramme N° 1 (comparer à photo 2). 


Oscillogramme n° 2 (comparer à figure 2). Remarquer la ‘‘distorsion’” 
sur la première demi-période. 

Oscillogrammes relevés au laboratoire (boucle de même caractéristi- 
que que celle proposée dans le programme). 
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Oscillogramme n° 3 (instabilité de fonctionnement autour de la fré- 
quence de capture). 


2FL ! 
fo 


Courbe de transfert relevée de la boucle (V2 =f{f;)) 
Oscillogramme n° 4. 





Annexe 1 


CHOIX DE VALEURS 
NORMALISEES 
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RECHERCHE DE VALEURS NORMALISÉES 


Lorque l'on a à trouver la valeur de composants électroniques obéis- 
sant à une certaine valeur, cela est très facile. Dans le cas d’une for- 
mule du type : 


ni 


RC 


(cas du décodeur 567), connaissant f, il est aisé, en se fixant R, d'en 
déduire une valeur de C. En général (loi dite de Murphy), cette valeur 
de C n'est pas normalisée. On se fixe alors une nouvelle valeur de R, 
d'où l’on déduit une nouvelle valeur de C, qui, elle aussi, n’est pas nor- 
malisée… Le problème est donc, ayant choisi une valeur normalisée pour 
R, d'en trouver une pour C, ou éventuellement la totalité des couples 
R, C normalisés obéissant à la relation. 

La capacité de calcul d'un micro-ordinateur est ici la bienvenue. Il est 
possible d'examiner par programme toutes les valeurs normalisées de 
R, puis, pour chacune d'entre elles, s'il existe, avec une certaine préci- 
sion, une valeur normalisée correspondante pour C. 

Nous supposerons que R peut prendre toutes les valeurs de la série E 
12, comprise entre 100 { et 1 MQ.C sera également pris dans la série 
E 12 et pourra prendre des valeurs comprises entre 1 nF et 1000 uF. 
Ces contraintes seront facilement modifiables. 


LES PROGRAMMES 
Ils ont tous la même architecture. 


Programme n° 1 


Lignes 50 et 60  :initialisation des valeurs normalisées. 
Lignes 110 à 150 : calcul de toutes les valeurs possibles de résistan- 
ces. 


Lignes 160 et 170 : pour chaque valeur de résistance, on cherche la 
valeur de C correspondant à la relation : 


_ 1 


RC 


Les valeurs trop petites (< 1 nF) ou trop grandes 
> 1000 nF) sont éliminées. 

Lignes 180 à 200 : pour une valeur de C convenable, on cherche la 
puissance de 10 par laquelle il faut multiplier C 
pour trouver un nombre supérieur à 1 (donc 
compris entre 1 et 10). 

Lignes 210 et 220: ce noveau nombre est comparé à +5 % aux dif- 
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férents nombres composant la série dans laquelle 
est pris le condensateur C. 

Lignes 240 et 250: s'il existe une valeur convenable, celle-ci est affi- 
chée. 

Ligne 270 : data contenant les 13 nombres de la série E 12 
dans laquelle sont supposés pris condensateurs et 
résistances. 


Le programme n° 2 distingue la série E 12 choisie pour les résistances 
(data en ligne 350) de la série E 6 choisie pour les capacités (data en 
ligne 360). 


Le programme n° 3 permet de calculer les éléments d’un multivibra- 
teur astable réalisé à l’aide d’un Timer type 555. La formule utilisée 
est alors la suivante : 


1,44 
(RA +2 RB)C 


Il'est demandé à l'utilisateur la valeur de f choisie et une valeur possi- 
ble pour RA. Le programme en déduit RB et C. 


LES VARIABLES 


V(13) : tableau contenant la série (E 12) des résistances. 

W(13) : tableau contenant la série (E 12 ou E 6) des condensateurs. 

R(52) : tableau contenant toutes les valeurs possibles normalisées de 
résistances (de 100 9 à 1 MQ). 


N : valeur de la puissance de 10 que contient le nombre représen- 
tant la capacité du condensateur. 

K : variable de défilement de toutes les valeurs de résistances. 

A : variable de défilement des valeurs normalisées de condensa- 


teur. 
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PROGRAMME 1 


F= oh R et C dans la série E 12 


RC 
À REM RERO RRRRERRRRRRERRERRE 
Z REX & CALCUL Li 
& REM # * 
3 REM # DE VALEURS NORMALIGEES # 
& REM # * 
7 OREM HHHKARERERERRERERERERREREE 
10 CL 
15 PRINT:PAINT" Ce programme calcule les valeurs" 
29 PRINT:PRINT*normalisees des composants R et C*° 
30 FRINT:PRINT*corresogndant a la formule f=i/8c,* 
33 PRINT:PRINT'R et C etant pris dans la serie E£° 
40 DIM VIA) :DIM NEJ):DIM R(S2) 


50 FGR I=f TO 13 

60 READ U:V(T)=UiWEII=U:NEXT ! 

70 PRINT: INPUT ‘Freguence=":F 

80 CLS:PRINT* F= *:F;"H2 "2: PRINT 
106 K=1:P=1 

116 FOR N=Z TG 5 

120 FOR I=i TO 13 

130 RIKI=V(I) {CAN 

140 K=K+i 


150 REXT LENEXT N 
170 C={/{F#R(K)):D=C: IF CO.000000001 OR Ci0.001 THEN 260 


180 N=0 

190 D=D#10:N=N+1 

200 iF Dè=t THEN 210 ELSE 190 

210 FOR Az T0 13 

220 IF D<1.05#W(A) AND D30.95#4{A) THEN 240 

230 NEXT A:G0TO 260 

240 PRINT'R=" ERIK): OHMG" st Ca" sM(A) "#0": 6-Ni micro." 
260 NEXT E 

270 DATA 1,1.2,1.9,1.8,2.2,2.7,3.3,2,9,4,7,5.6,6.8,8. 2, 10 
280 G0TO 280 

8900 END 

9000 2=0 

9100 REM #COPIE ECRANX 

9110 RESTORE 9199: T0T=0 
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9140 MEMORY &SFFF 

9150 FÜR ï= &AGDO TO &AÜEO 
9160 READ V:POKE IV: TOT=TOT#V 
9170 NEXT I 


0180 IF TÔTES 20364 THEN PRINT “ERREUR EN DATA" ELSE CALL £&ACOOC:RETURN 
9190 DATA 205,{66,1 Le 62,27,205,157,182 

9200 DATA 62 ,49, 2 5,157, 160,205, 106,167 

9210 DATA : 205,231, 187, 50, 180, 60, 17,8 


220 DATA 0,33,143,1,34, 178, 160, 82 
225 DATA 7,50,177,160, 62, 10,205, 157 
9230 DATA 160,0,0,0,0, 0,62, 27 
9240 DATA 205,157, 160, 62, 76,205, 157,160 
9250 DATA 62,127,205,157, 160, 62, 2, 205 
9260 DATA 157,160, 14,0,58, 177, 160,71 
9270 DATA 229,197,213, 205,240, 187, 209,193 
9280 DATA 33,180, 160, 190, 225,55, 32, 1 
9290 DATA 167,203,17,43,43,16,233,58 
9300 DATA 177,160,254,7,40,7, 175,203 
9310 DATA 17,203, 17,203, 17,121,205,197 
9320 DATA 160,19, 229, 33, 127,2,55,227 
9330 DATA 82,225,56,5,42,178, 160,24 
9340 DATA 193,25, 124,181, 40,32, 43,17 
9350 DATA 0,0,36,178, 160,62, 7, 189 
9360 DATA 32,146,124, 180,32, 142,62,4 
9370 DATA 50,177, 160,24, 135,205, 46, 169 
9380 DATA 56,251, 205,43, 189, 201, 62,27 
9390 DATA 205,157,160,62,64,205, 137, 160, 201 


F= 11509 “1z 
RB= 129 OHMS L= S.6 “IDT-1 micro-r. 
R= LE CHMAS C= 4.7 ÆH1NT-1 micro—-s. 
R= 229 OHMS C= 5.9 #10f-1 micro-r. 
R= 2719 OHA£ C= 3.23 æ#19f-1 micro-r. 
H= SED Oh'MS C= 2.7 H#IAOT-1 micro-F. 
R= 229 OFM£S C= 2.2 #I19OT-1 micro-E. 
R= 470 OMS C= 1.85 #107-1 micro-r. 
KH= 560 OHMS C= 1.5 #101t-1 micro—-F. 
R= 15509 OHmS C= $S.5 #i0Ot-2 micro-!. 
B= 1569 CHAMS C= 34.7 #I1OT-2 micro—-r. 
R= e2@9 OMS C= 2.9 æ1G0f-2 micro-F. 
E= LV OHMS C= 3.3 #IOT-2 micro-F. 
P= 2297 CHMS C= 2.7 #æ1i0f-8 micro-xr. 
R= 3260 OMS C= 28.8 w#iDt-e micro-r. 
R= 4709 OHAS C= ÎL.Ë #iOT-EÉ nicro-r. 
R= S66 Ohms Cz 1.5 xibt-e micro-r. 
R= 15868 CHMS = S.5 #197-35 micro-r 
R= 156009 OHMS C= 4.7 #1QT—-Z micro-s 
R= LEO OHMS£ C= 2.92 107-232 micro-r 
= 27999 OHMS C= ©. #109T—-2 micro-r 
R= Z22:298 CHMC C= 2.7 HICT-S nicro-r. 
Ez= 355699 OCrMS £= 2.2 K1iSGt-Z5 nicro-r. 
Az 47999 OCHMS C= 1.3 #iO0Of-2 micro—-z. 
H= 5956099 CHMS C= 1.5 Kigÿt-z nmicro-r. 
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Fâäz £509 ORMS C 4.7 #18T-2 micro-F. 
Rf= 18560 CHMS C 3.3 #1O0T-Z micro-F. 
R= 393699 OHMS C= 2.2 #1Ot-3 micro-Fr. 
RR= 106606 OHMS C pu 
R$= 180005 OHMS C 1 


Æ1OT-T3 micro-r. 


K1HT—-5 micro— 


F=  £2880 Hz 
RB= 5500  OKMS 

PA= i5nn OHMS C= 28.2 #i6t-3 micro-F. 

RA= 12569 OHMS C= 2.8 #i0t-2 micro-F. 

RA= 28208 ORME ©= 2.5 #iGt-5 micro-F. 

Rfiz Set ONMAS CL 1.S K1GT-5 micro—-Fr. 

1C-REM HH#RRARRRREREERÉERERSRESS 

2G REM + CALCUL * 

J0 REM # 4 

40 REM + DE VALEURS NORMALISGEES # 

SÛ REM * * 

60 REM #RRI RERO ER ÉREE 

70 CLS 

BÛ PRINT:PRINT" Ce programme caicule les valeurs 


90 PRINT:PRINT'rormalisees des composants À et C° 

100 PRINT:PRINT" correspondent a la formuie f=!/R€, 

110 PRINT:PRINT'R etant pris dans la serie FiZ et © dansia seri 
120 DIM V(19):DIM W(1G):DIM R(52) 

130 FOR I={ Tü IJ:READ YJ:VITIU:NEXT I 

140 FÜR I=1 T0 7:READ U:HII)=U:NEXT : 


C 3 
re 


A 
Fa 
Le 


" 





O PRINT:INPUT "Frequence=":;F 
160 CLS:PRINT" F= “F3; *Hzt: PRINT 
180 K=1:P=i 
190 FOR N=2 TO 5 
17 


200 FOR I=i TG 13 
210 RIK)=V(TI#IDAN 
220 K=K+#i 
230 NEÏT 
240 FOR K=1 79 53 

230 C=1/ CFRRCKS D:O=C:IF CéG.0G0G000Gi GA CC. CC1 THEN 330 
266 N=G 

2170 D=D#IO:N=N+! 
280 IF Di=i THEN 290 ELSE 279 
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NEXT N 
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R A=1 TO 7 


306 IF Dii.05#ki 
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350 DATA 1,1. 
360 DATA 1,1.5 
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10 REM HER RER ER END EERE 


20 REM # CALCUL # 
30 REM + 4 
40 REM +# L 
50 REM *# DE VALEURS NÜRMALISEES # 
60 REM # Li 


70 REM KFAREFRARESRRREEFÉFÉREEREES 
Bo LS 


90 PRINT:ORINT" Ce programme calcule les valeurs* 
100 PRINT:PRINT'ncrmalisees _ composants KR et £° 
110 PRINT:PRINT'correspordant à la formule": PRINT 
BjC.® 

115 PRINT:PRINT'RA etant pris dans la serie Elite 


Eb,* 


117 PRINT:PRINT'F et RE etant donnees, 
120 DIM V{ig): Di W(1S):DIM R(52) 
130 FOR I=1 TG 13:READ U:VIT)=iÜsNE) 
140 FOR I={ TO 7:READ U:WiTI=U:NEX 


150 PRINT: AE “Frequerce=":f 
155 PRINT:INPUT "Valeur de RE:":FE 
180 CLS:PRINT" 

185 PRINT* 

180 K={:P={ 

190 FÜR N=Z Tü 5 
206 FOR I=1 TG 13 
219 RiK)=V{I)#IOMN 
220 K=K+i 

230 NEXT IiNCXT N 
240 FOR K=1 Tü 57 
250 C=1,44/(F#{RIKI4CARB) 3: DCS IF 
Z50 N=0 

270 D=L#{0:N=N+{ 

280 IF Li=1 THEN 290 ELSE 270 

290 FÜR Ai TO 7 

360 IF D£f.0S#4{A) AND D:0.954W{A 
Ji0 NEXT A:GüTC 349 
320 PRINT:PRINT®RA= 
iCA*s6-Nitmicro-F.* 


FiRCE) TAR: 


3 NEXT k 

350 GÈTÉ 35 
S60 DATA 1,1.2,1.5,1,8,2,2,2,7,3, 
370 DATA 1,1,5,2,2,4,9,4,7,6,8,i0 
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F= NU # 
RE= "iRBst GENS": 
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PRINT f= 


4 
à 


THEN 349 


TAF (20) 
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lu 
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a 
[HI 
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ol 
3 
3 
it 


ki 


& 
composants K et © 
e 








normaiisees dez 
sorrespondant à 1x formul 
f=i.A14/ (RA+ENREDIC. 
EN se 1x serie Elec et C 
F eat RE etant donneszs. 
e=? 115939 
e RB:7 16565 
F= 11509 Hz 
RB= 168955 CGARMS 
R$= SsS59 GHMS C= 4,7 #iNt-3 micro-F. 
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